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CHAPTER I. 



INTRODUCTORY, DEFINITION OF SPHERICAL HARMONICS, 

1. If V be the potential of an attracting mass, at any 
point Xy y, z, not forming a part of the mass itself, it is 
known that Fmust satisfy the differential equation 

d^V d^V d^V ^ 

^+3^ + ^ = ^ W' 

or, as we shall write it for shortness, V*F=0. 

The general solution of this equation cannot be obtained 
in finite terms. We can, however, determine an expression 
which we shall call F,, an homogeneous function of a;, y, z 
of the degree i, i being any positive integer, which will 
satisfy the equation ; and we may prove that to every such 
solution Fj there corresponds another, of the degree — {i + 1), 

F 

expressed by -^^ , where r* = a?' + y* + z\ 

For the equation (1) when transformed to polar co-ordi- 
nates by writing a; = r sin ^ cos ^, y = r sin ^ sin (j)yZ = r cos 0, 
becomes 



d'irV)^ 1 df. .rfF\ 1 cPF ^ ,^,. 

And since V satisfies this equation, and is an homo- 
geneous function of the degree i, V^ must satisfy the equa- 
tion 

F. H. I 
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siDce this is the form which equation (2) assumes when F 
is an homogeneous function of the degree /. 

Now, put F] = r**** Z7], and this becomes 

or 
i(*>l)tr, + -i-.|.fsin^§UJ^5? = (2). 

Now, since 27] is a homogeneous function of the degree 

- (i + 1). 

dr * dr 

= -iU,', 
dr' dr 



or 



r^P^m+DU,: 



therefore equation (2) becomes 



d'(rU) 
r- — — i- + 

ar" sin dff 



1 d f . .dU\^_J_drU,_ 



shewing that U^ is an admissible value of F, as satisfying 
equation (2). 

It appears therefore that every form of Z^ can be ob- 
tained from F,, by dividing by r'**\ and conversely, that 
every form of F, can be obtained from U^ by multiplying 
by r*'*\ Such an expression as F, we shall call a Solid 
Spherical Harmonic of the degree i. The result obtained 
by dividing F, by >*, which will be a function of two inde- 
pendent variables and ^ only, we shall call a Surface 
Spherical Harmonic of the same degree. A very important 
form of spherical harmonics is that which is independent 
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of <j>. The soKd harmonics of this form will involve two of 
the variables, x and y, only in the form a? + i^, or will be 
functions of a;* + 3/* and z. Harmonics independent of <f> are 
called Zonal Harmonics, and are distinguished, like spherical 
harmonics in general, into Solid and Surface Harmonics. 
The investigation of their properties will be the subject of 
the following chapter. 

The name of Spherical Harmonics was first applied to 
ihese functions by Sir W. Thomson and Professor Tait, in 
their Treatise on Natural Philosophy. The name " Laplace's 
Coefficients" was employed by Whewell, on account of 
Laplace having discussed their properties, and employed 
them largely in the Mecanique CMeste. Pratt, in his 
Treatise on the Figure of the Earth, limits the name of 
Laplace's Coefficients to Zonal Harmonics, and designates 
all other spherical harmonics by the name of Laplace's 
Functions, The Zonal Harmonic in th6 case which we shall 
consider in the following chapter, i.e., in which the system 
is symmetrical about the line from which is measured, 
was really, however, first introduced by Legendre, although 
the properties of spherical harmonics in general were first 
discussed by Laplace; and Mr Todhunter, in his Treatise, 
on this account calls them by the name of "Legendre's 
Coefficients," applying the name of " Laplace's Coefficients " 
to the form which the Zonal Harmonic assumes when in 
place of cos 6, we write cos 6 cos ff + sin 6 sin 6' cos (^ — ^'). 
The name " Kugelfunctionen " is employed by Heine, 
in his standard treatise on these functions, to designate 
Spherical Harmonics in general. 



\— *! 



CHAPTER 11. 



ZONAL HARMONICS, 



1, We shall in this chapter regard a Zonal Solid Har- 
monic, of the degree i, as a homogeneous function of 
(«* +y^^f and z, of the degree t, which satisfies the equation 

d^^ dy""^ dz^^ ' 

Now, if this be transformed to polar co-ordinates, by 
writing r sin 6 cos ff> for a;, r sin 6 sin <^ for y, r cos for z, we 
observe, in the first place, that a;'+y" = r*sin*5. Hence 
V will be independent of <j>, or will be a function of r 
and only. The differential equation between r and 
which it must therefore satisfy will be 

ar sm5 rfd V d0J 

Now y, being a function of r of the degree i, may be 
expressed in the form r^Pi, where P, is a function of only. 
Hence this equation become^ 

or, putting cos0=fi, 

In accordance with our definition of spherical surface 
harmonics, P, will be the zonal surface harmonic of the 
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degree l When it is necessaiy to particularise the variable 
involved in it, we shall write it P^ (ji). 

The line from which is measured, or in other words 
for which /a = 1, is called the Axis of the system of Zonal 
Harmonics; and the point in which the positive direction 
of the axis meets a sphere whose centre is the origin of 
co-ordinates, and radius unity, is called the Pole of the 
system. 

Any constant multiple of a zonal harmonic (solid or 
surface) is itself a zonal harmonic of the same order. 

2. The zonal harmonic of the degree i, of which the 
line /A = 1 is the axis, is a perfectly determinate function of 
fi, having nothing arbitrary but this constant. For the 
expression r^P^ may be expressed as a rational integral 
homogeneous function of r and ^, and therefore P^ will be 
a rational integral function of cos 0, that is of fi, of the 
degree i, and will involve none but positive integral powers 

But P^ is a particular integral of the equation 



d/M 



and the most general form oi fQi) must involve two ar- 
bitrary constants. Suppose then that the most general 

form o{f(ji) is represented by P^ jvd/M. We then have 



+ 2(l-^«)g. + P.|{(l-^«)..}. 



Hence, adding these two equations together, and ob- 
serving that, since P, satisfies the equation (3), the coefficient 
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of Ivdfi will be identically eq[ual to 0, we obtain, for the de- 
termination of V, the equation 

whence i',a-/**)^+2{(l-/**) ^-/*i*j« = 0, 



or 



iFHif-i^-)^"-"- 



the integral of which is 

log V + log P* (1 — ytt*) = log Ci = a constant ; 

_ (7, 

Hence jvdfi ^^C+C.j p, n _»,''; '* 

and we obtain, for the most general form oif(ji), 

Now, P, being a rational integral function of fi of i 

dimensions, it may be seen that Itz ^rrm will assume the 

form of the sum of i + 2 logarithms and i fractions, and 
therefore cannot be expressed as a rational integral function 

of /Aw Expressions of the form P^ Ij- ,^ p^ are called KugeU 

functionen der zweiter Art by Heine, who has investigated 
their properties at great length. They have, as will hereafter 
be seen, interesting applications to the attraction of a sphe- 
roid on an external point. We shall discuss their properties 
more fully hereafter. 

3. We have thus shewn that the most general solution 
of equation (2) of the form of a rational integral function of a 
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involves but one arbitrary constant, and that as a factor. 
We shall henceforth denote by P„ or P^ (fi), that particular 
form of the integral which assumes the value unity when fj, 
is put equal to unity. 

We shall next prove the following important proposition. 

Ifh.be less than unity, and if {I- 2/Lth + h*)"*^ he expanded 
in a series proceeding hy ascending powers of h, the coefficient 
ofhf will ftePp 

Or, (l«2/ih + hT*J = Po + P|l^ + - + Pili*+..- 

We shall prove this by shewing that, if H be written for 
(1 — 2fJb + A*)"*, H will satisfy the differential equation 

For, since E=> (1 - 2M + A')-* 

h dfi ' 



= - 2/tff' + 3 (1 - fi') H' ^^ 

«5' {1 - 2M +^' + A (m- A)} 
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« 

= 2 (I- fill) E'(ji-h) -fiH*. 

= ^ 3/.ir' + 8 {(1 - A**) A + (1 - M) O - ^)} -S" 
= -S,iE»+3[/i{l+h')-2/M'h}H* 

= - 3/* {S' - (1 - 2M + h*) H'} 
= 0, since 1 - 2^A + ;i'' = ^^. 

Therefore |{(1 - ^^ ^} + aJ[,(AF) = 0. 

This may also be shewn as follows. 

If a?, y, z be the co-ordinates of any point, / the distance 
of a fixed point, situated on the axis of i, from the origin, 
and R be the distance between these points, we know that. 



and that 



is-»- 



Now, transform these expressions to polar co-ordinates, 
by writing 

a; = rsin^cos^, y = rsin^sin^, j? = rcos^, 

and we get 

and the differential equation becomes 



dr 



''©+ii^^h^^©}=^' 
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or, putting cos 5 = /i, 

Now, putting r=z'h, we see that 



= 0. 



or 



IH 
B~ z" 

T 



and 



dAw^z'dh*^^^' 



.'. the above equation becomes 



or 






4. Having established this proposition, we may proceed 
as follows: 

If p^ be the coefficient of h* in the expansion of -ff, 
11=^1 ■^pJb+pJi!' + ...-\-pJi* + ... 



• • • 



(P 



Also, the coefficient of h* in the expansion of 

Hence equating to zero the coefficient of h\ 
|{C-''')|'}+''(.-H-l)?.-«. 
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Also Pi is a rational integral function of fi. 
And, when /t* = 1, iZ'= (1 - 2A + h^^ 

Or when /a = 1, jp, = 1. 
Therefore p^ is what we have already denoted by P,. 
"We have thus shewn that, if h be less than 1, 

If h be greater than 1, this series becomes divergent. 
But we may write 



''hV''^f'^''''^h''^"')' 



since , is less than 1, 
h 



P P P 

£?4.-V+ 4.-1* 



-s 1-2 4- - Ji 4- 4. -_* 4- 



Hence P, \s also the coefficient of A-^^+^> in the expan- 
sion of (1 — 2fih + /t*) "^ in ascending powers of r when h is 

greater than 1. We may express this in a notation which is 
strictly continuous, by saying that 

This might have been anticipated, from the fact that the 
fundamental differential equation for P, is unaltered if 
— (i + 1) b6 written in place of i ; for the only way in 
which % appears in that equation is in the coefficient of P„ 
which is % (i + 1). Writing — (t + 1) in place of t, this be- 
comes — (i + 1) {— (f + 1) + 1} or (i + 1) i, and is therefore 
unaltered. 



5. We shall next prove that 

where 7^ = 0? + ^-+ z\ 

Let l={^ + y' + z*)-^'=f{z), 

and let k be any quantity less than r. 

Then [x^ + y^ + (z-' ky}'^ =/ (« - k), 

and, developing by Taylor's Theorem, the coefficient of k' is 

(.iyJ^..or(-l)'-J— .|.(1). 
^ '1.2... I ^ ' 1.2 ...I a^*\.r/ 

Also {a? + 2/* + (« - A;)»ri = (r» - 2i!;« + A;')"* 



1/, - & ifc»\-* 



Since z = /tr, 

in the expansion of which, the coefficient of k* is 

, p. 

Equating these results, we get 






The value of P, might be calculated, either by expanding 
(1 — 2fjbh + h^"^ by the Binomial Theorem, or by effecting the 

differentiations in the expression (— l)'n— ir~o •Tnf-jj 

and in the result putting - = /t. Both these methods how- 
ever would be somewhat laborious ; we proceed therefore to 
investigate more convenient expressions. 
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6. The first process shews, by the aid of Lagrange's 
Theorem, that 

1 /7* 

Let y denote a quantity, such that 

h being less than 1. 

Then 

1 
c/y h 



dfi 



/l_2^+l^* (1-2M + A')4' 



Also (y-|J = l-^ + p; 

*• ^ A -^ A ' 
Hence, by Lagrange's Theorem, 

therefore, differentiating with respect to n and obserying 
that 

|=(l-2M + A')-i. 
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7. From this form of P^ it may be readily shewn that 
the values of fi, which satisfy the equation P^ == 0, are all real, 
and all lie between — 1 and 1, 

For the equation 

(/A*— 1)* = has 4 roots =1, and i roots = — 1, 

f\ -p {fjf-- 1)*= has f — 1 roots = 1, (t — 1) roots = — 1, and 
one root = 0, 

-T-5(/Lt* — 1)*«0 has (t — 2) roots = 1, one root between 1 

and 0, one between and = — 1, and (i — 2) roots = —1, 
and so on. Hence it follows that 

d* i i 

i-j (jjf — 1)* = has ^ roots between 1 and 0, and - roots be- 
tween and — 1, if % be even, 

i— 1 i— 1 

and —^— roots between 1 and 0, — ^— roots between and 

1, and one root = 0, if t be odd. 

It is hardly necessary to observe that the positive roots of 
each of these equations are severally equal in absolute mag- 
nitude to the negative roots, 

8. We may take this opportunity of introducing an im- 
portant theorem, due to Rodrigues, properly belonging to 
the Differential Calculus, but which is of great use in this 
subject. 

The theorem in question is as follows: 
Ifvsa he any integer less than i, 
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It may be proved in the following manner. 

» • 

If {a? — 1)* be differentiated t - m times, then, since the 
equation 

(x*-l)' = 

has i roots each equal to 1, and % roots each equal = — 1, it 
follows that the equation 



djo 



has i — (i — m) roots (i. e. 7n) roots each = 1, and m roots 
each = — 1, in other words that {a? — 1;"* is a factor of 



We proceed to calculate the other factor. 

For this purpose consider the expression 

(a? + aj (a?+ aj ... (aj + a,) (a; + i8J (aj + /3^ ... (a? + A). 

Conceive this differentiated (I) i — m times, (II) i + m 
times. The two expressions thus obtained will consist of an 
equal number of terms, and to any term in (I) will corre- 
spond one term in (II), such that their product will be 
{x + a^ {(v + a^) ... (a?+aj {x + ^^) {x + ^^) ... (a: + i8,),i.e. the 
term m (II) is the product of all the factors omitted from 
the corresponding term in (I) and of those factors only. 
Two such terms may be said to be complementary to each 
other. 

Now, conceive a term in (II) the product of p factors of 
the form x + a, say a;+a', x + a ...cc + a^', and of ^factors 
of the form a;+^, say a; + ^^, aj + i8^^...aj + ^jjj. We must 
have J!? + J = i — m. 

The complementary term in (I) will involve 

p factors x + ^, x + fi" ...x+^\ 

q factors x + a^, a? + a^, . . . x + Oj^,. 
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Now, every term in (I) is of i + m dimensions. We have 
accounted for p + j (or i — m) factors in the particular term 
we are considering. There remain therefore 2m factors to 
be accounted for. None of the letters 

a, a ...a'', p^> p^, ...P(j), 

can appear there. Hence the remaining factor must involve 
VI a's and m /3's, — say, 

There will be another term in (II) containing 

(0, + ^) {x + ^") ... {x + ^') {x + a) (x + a,) ... (a? + a,J. 

The corresponding term in (1) will be, as shewn above, 
(x + a:){x + a:')...{x + a!^){x + l3){x + ^,:)...{x + ^J 
(x + ^a)(x + ,^)...{x + „i){x + ,fi){x-h^)...{x + ^). 

Hence, the sum of these two terms of (I) divided by the 
sum of the complementary two terms of (II) is 

(x + ,fi) (x + ^a) ... (x + ^oi) {x + ,&) (ic + ,/3) ... (a;+j9). 

Now, let each of the a's be equal to 1, and each of the ffs 
efjual to — 1, then this becomes (a^ — l)'^ The same factor 
enters into every such pair of the terms of (I). Hence 

Or ^^,-^ = (ai» - 1)*" ^1^^ , to a numerical 

factor prh. 

The factor may easily be calculated, by considering that 
the coefficient of x*"^ in X*"^ is2i(2i-l)...(i+m+l), 

and that the coefficient of x*"^ in zri+^i — '^^ 

2i (2i- 1) ... (i + m+ 1) (i + m)... {i-m'^' 1). ' 
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Hence the factor is 

1 1 . 2 . . . (i — m) 

{i+m)[i+m-l)...{t-m + iy^^ 1.2... {i + m)' 

9. This theorem affords a direct proof that (7-7-j (/a*— 1)', 

C being any constant, is a value of / (fi) which satisfies the 
equation 



i (('-"•) 



^)}+.-(.-+i)/W=o. 



=.-(.-+l)^.0.'-l)' 

from above, 
or 

I [(>-"•) ^f|. t'- - ')')] +<(.>1){|0'--V}=o. 

Hence, the given differential equation is satisfied by put- 

Introducing the condition that P^ is that value of / Qi) 
which is equal to 1, when /* = 1, we get 

•^* " '2^1723- 5;? ^*"'^^'' 

10. We shall now establish two very important proper- 
ties of the function P^ ; and apply them to obtain the develop- 
ment of P« in a series* 
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The properties in question are as follows : 
Ifi and m be unequal positive integers, 

And Ij'^^-^V 

The following is a proof of the first property. 
We have 

Multiplying the first of these equations by P^ the second 
by Pp subtracting and integrating, we get 

+ {i (i+ 1) - m (m + 1)} jp.P^dfi = 0. 

Hence, transforming the first two integrals by integration 
by parts, and remarking that 

i{i+l) — m(m-f 1) = (i — m) (i+m + 1), 
we get 

+ (i -m){i + m + 1) fPtPJii = 0, 



or 



(1 -,*«) (p„^'_P,^") + (i-m)(t + m + l)/p^,<f^=0, 

since the second term vanishes identically. 

P. H. 2 
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/: 



Hence, takijbg the integral between the limits — 1 an4 
+ 1, we remark that the factor 1 —/a* vanishes at both limits^ 
and therefore, excqpt when i — m, or « + m + 1 = 0, 

-1 
We may remark, also that we have in general 

a result which will be useful hereafter. 

11. We will now consider the cases in which 

i — m, or i + m + 1 = 0. 

We see that i-f m + 1 cannot be equal to 0, if i and m are 
both positive integers. Hence we need only discuss the 
case in which m = i. We may remark, however, that sjnc^ 

Pi = P_(i+i), the determination of the value of I PJ dfi wiU also 

•^ -1 

give the value of I P^P^^^^)dfA. 

The value of I P^d/M may be calculated as follows : 

{l-^2fih'hh^)'^ = P, + PJi + ,.. + PJi*+...; 

.-. {1^2H^h + h')-' = {P,+ PJi+... +p,h' + ...y 

= Po' + P,W+...+P,T'+... 
+ 2P,PJi + 2P,PJi' +... + 2P,PJi' + . . . 

Integrate both sides with respect to fi ; then since 

j{1^2fih + hYdfi^'-^log(1^2fih + h'), 
we get, taking this integral between the limits — 1 and + 1, 

Jl0g^ = J' P.'dfl + h'j"^ P,'dfJL+...+h''j'^ P,V/*+... 

all the other terms vanishing, by the theorem just pio/ed. 
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Hence 2fl + g- + ... + q- . t +>" J 
Hence, equating coefficients of A**, 



/ 






12. From the equation J PJP^dfi = 0, combined with 

the fact that, when /a = 1, P^ = 1, and that P^ is a rational 
integral function of /a, of the degree i, P, may be expressed 
in a series by the following method. 

We may observe in the first place that, if m be any 

integer less than i, I yU^P^dii = 0. 

For as P^, P^.!... may all be expressed as rational in- 
tegral functions of /a, of the degrees m, w — 1 . . . respectively, 
it follows that yS^ will be a linear function of P^ and zonal 
harmonics of lower orders, /a"*"* of P^j and zonal harmonics of 

lower orders, and so on. Hence l/x*"P,(i/A will be the sum of 

a series of multiples of quantities of the form I P„P<c//a, 

m being less than i, and therefore I yS^^Pfi^L = 0, if m be any 
integer less than i. 
Again, since 

(l-2/ife + AV* = Po + P,A+...+P,A,'+... 
it follows, writing — \ for A, that 

(i + 2M+Ar*=-Po-P»A+.,. + (-iyp,/t'+... 

2— a 



20 ZONAL HABMONIC?S. 

And writing — /a for /a in the first equation, 

P^', P^,.,PI ... denoting the values which P^, P^...Pi, 
respectively assume, when —/a is written for /a. Hence 
PI=iP^ or —Pi, according as i is even or odd. That is, 
P^ involves only odd, or only even, powers of i, according 
as i is odd or even*. 

Assume then 

Our object is to determine A^y A^_^ .... 

Then, multiplying successively by /tt*~*, /a*"*, ... and inte- 
grating from — 1 to + 1, we obtain the following system of 
equations : 



!- 4- 4- «-» 4. s= 



2i-l 2i-3 2i-2s-l 



"f" 77": ^ "T" • • • "T ^ « T^j ^ T" • • • "~ V, 



2i-3 ' 2i-5 2i-2«-3 



-^< ■ I -^<-« I I -^<-2f 



+ o» O o + '**+^ P^ + ... = 0. 



2i_25-l ' 2i~25-3 2i-4<s-l 



And lastly, since P, = 1, when /* = 1, 

•^i "T" -"-i-j "T • • • T -"i^ji T • • • ^^ -I- 5 

the last terms of the first members of these several equa- 
tions being 

AAA 

-r— V* -;— ^...-T^, A, if ibe even, 
t — i t — o i 

A A A 

-7— *7j, -; — V-**-^* ^i» if ^be odd. 
t — 2 i — 4 2 ^ 

13. The mode of solving the class of systems of equa- 
tions to which this system belongs will be best seen by 
considering a particular example. 

* This is also eyident, from the fact that P« is a constant multiple of 
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Suppose then that we have 



a + a 6 + ot q+ a 
X y z 



= 0, 



= 0, 



X y , z 1 



a + G) 6 + fti c + o) o) ' 

From this system of equations we deduce the following, 
6 being any quantity whatever, 

X y ^ ^1 (g-g)(g~i8)(a+a>)(& + a>)(c+a)) 

a+6'J + 6^"^c + ^ ft) (fti-a)(a)-/3)(a+5^(6+^)(c+5)" 

For this expression is of — 1 dimension in a, h, c, a, /9, 7, 
6, ©; it vanishes when ^ = a, or = fi, and for no other 

finite value of 0, and it becomes = — , when = 0). 

We hence obtain 

Tx(.xfi\( y I M_l (g-Q^)(g-/3)(«+^)(&+^)(c+a) ) 
^"^r"*"^^V6+^ c+6^y "" o> (a,-a)(fl>-/3) (6+6^)(c+6;) ' 

and therefore, putting ^ = — a,^ 

1 (a + a)(a + /9) (a + co) (6 + a)) (c + g)) 
ft) (a — ^^a — c) (ft> — a) (o)— )S) ' 

with similar values for y and «, 

And, if ft> be infinitely great, in which case the last 
equation assumes the form a? + y + ^r = 1, we have 

_ {a + a) (a-hff) 
^^ (a-6)(a-.c)' 

with similar values for y and z. 

14. Now consider the general system 

-5-+-^^+...+-^ + ...=0, 

*' .— ^'-. +...+ ^'7* +... = 0, 



tti + a*^ »w+a«-» ««-«. + «*-t 
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"^^ - '^'-- + ...+ — ^^=^^— +...=0, 



tt< + ««-8f C^i-a + ««-a» tt«-» + «<-M 



""' +^^ + ...+-^^=!^ + ... = -; 



a, + G) a,_,+ o> ^i^ + co ft) 

the number of equations, and therefore of letters of the 

forms X and a, being — ^ if i be odd, ^ + 1 if i be even ; and 

i — 1 . 

the number of letters of the form a being — ^— if t be odd, 

and ^ — 1 if t be even. 
We obtaio, as before. 



^i-W"~ 



_ 1_ (g-g< ) (^-0>»»(^-gi-J >*' (g, + o>) (a,,,+ ft)) . > . (a<,a.-fft)) . . . 

and, multiplying by a^^^ + ^, and then putting tf = — a,_^, 

ft> (fi)^aj (to-at,^ ... (ft) - a,.J ... 

(tt j + ft)) (a<_,+ ft))...(a,^a,+ft))... 

15. To apply this to the case of zonal harmonics, we see, 
by comparing the equations for x with the equations for A, 
that we must suppose » = x ; and 

<^i = h a|_, = i — 2,...a,_g, = i- 2«... 

a, = t — 1, a|_j, = i — 3,...a|.j, = i — 2«— 1... 

Hence 

. _ (2z-2g-l)(2^-2g-3)...{2(t-2g)-l|.., 
*-*• (-2«) {-(25-2)}...{(i-2«-l) or (i-2d)} 

_. ., (2t-2g-l)(2i-2g-8)...{2(t-28)--l|... 
^ ^ aj(2«-2),..2x2.4.,.ii-r2«-l)or(i-25)' 
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Or, generally, if i be odd, 

' 2.4...(i-l) ' 

. (2t-3)(2f-5)...t 

'-»" 2.4...(i-3)x2 ' 

J _ (2t-5)(2t-7)...(z-2) 
'-*~ 2.4...(i-5)x2.4 ' 



, _(_-i\T t(t — 2)...3 
»~^ ^ 2.4...(i-l)* 

And, if I be even, 

^_ (2i-l)(2t-3)...(i+l) 
'^ 2.4...i 

. (2t-3)(2t-5)...(t-l) 

'-• 2.4...{i-2)x2 ' 

. _ (2i-5)(2i-7)...(t-3) 
*-* 2.4...(t-4)x2.4 ' 



We give the values of the several zonal harmonics, from 
P, to Pjo inclusive, calculated by this formula, 

■^» 2^^ 2' 



_ 3/t''-l 

~ 2 

5 , 3 

2'*~2 

5fi' - 3/[t 

2 



•P.^o/^'-o^ 



• « 
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p__7.o^ 5. O3O.X 

35/**- 30/** + 3 
8 ' 

» 2.4^* 2.2'* ■'■2.4!'* 

_ 63/t' - 70/t' + 15/1. 
~ 8 ' 

* '^ 2.4.6 '* ~ 2.4x:2'* ■•" 2 x 2 . 4^* ~ 270 

_ 23 V- 315/t*+ lOo/t*- 5 

16' 

„ 13.11.9 , 11.9.7 , 9.7.5 , 7.5.3 
' 2.4.6 '* 2.4 x2'* ■'■2x2.4'* 2.4.6^* 

_ 429/^' - 693/1" + 315/t' - 35/t 
15 • 

15.13.11.9 , 13.11.9.7 , 11.9.7.5 « 
*~ 2.4.6.8 '* 2.4.6x2'* ■'■2.4x2.4'* 

9.7.5.3 , ^ 7.5.3.1 
2x2.4.6'* ■ 2.4.6.8 

_ 6435/t' - 12012/t' + 6930/t* - 1260/t* + 35 

128 

P _ 17.15.13.11 , 15.13.11.9 T , 13^11j9^ . 
' ~ 2.4.6.8 '* 2.4.6x2 '*"'" 2.4x2.4 '* 

11.9.7.5 , 9.7.5.3 
2x2.4.6'* ■'■2.4.6.8'* 

_ 12155/t'- 25740/t' + 18018/t'- 4620/t' + 315/t 

128 ' 

p ^ 19.17.15.13.11 I, 17.15.13.11. 9 ,1 5.13.11.9.7 , 
" 2.4.6.8.10 '^ 2.4.6.8x2"'*"''2.4.6x2.4'* 
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18.11.9.7.6 ^ 11.9.7.5.3 ,_ 9.7.5.3.1 



2.4x2.4.6^ '2x2.4.6.8'^ 2.4.6.8.10 

_ 46189/i*"- 109395^^+ 90090/^"- 30030/1,*+ 3465/i^'-- 63 

256 

' It will be observed that, when these fractions are reduced 
to their lowest terms, the denominators are in all cases 
powers of 2, the other factors being cancelled by correspond- 
ing factors in the numerator. The power of 2, in the 
denominator of P„ is that which enters as a factor into the 
continued product 1 . 2...i. 

16. We have seen that / /i*" P, . e?/i = 0, if m be any 
integer less than L 

It will easily be seen that if w + i be an odd number, the 
values of I/a"' P^ . dfi are the same, whether /a be put = 1 or 

— 1 ; but if m + i be an even number, the values of Ifi"^ P^. dfi 

corresponding to these limits are equal and opposite. Hence, 
(m + i being even) 

then f /i^P^.d'/i^O, if m = «-2, ;-4 

Jo 



and 



We may proceed to investigate the value of I jjJ^P^. dfi, 

Jo 
if m have any other value. For this purpose, resuming the 

notation of the equations of Art. 13, we see that, putting 

tf s= m -f 1, and © = x , we have 

^J -j. ^*-3 f. 4- ?^ + 



a^+m + l a^^ + m + l a^^ + m+1 

_ (m + 1 — 0L^{m + 1 — g^) ... (m + 1 — o^_^^ ... ^ 
(a, + w+l)(a,^ + ^+l) ...(a,^ + w+ 1)...' 
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and therefore, putting a?, = J., . , ., a, = « . . ., a, = / — 1 . . ., 

we get 

T • • • "i *^S S ; rT T .• • • 



ri A 



_ (m-t+2)(m-t + 4)...(m-l) ^-i,,qjj 
""(w + i + l)(m + i-l)...(w + 4)(m + 2)'^*^®^''' 

, (m — 1 + 2) (m — 1*+4) .., m .« ., 

^d =/ — . > . Tx/ — . . Tx / . QN / — -^rr-if i beeven, 

(wiH- » + 1) (^n + i — 1) ... (m + 3) (m + 1) 

In the particular case in which wi = t, we get 
/, f'^*^'' = C2e + 1) (L-- 1) . . (^-^ 4) (e> 2) ^^'^'^'^^^ 

'^^ = (2/+l)(2t-i).'.'.'(>-+3)(i + l) (* ^^^^^- 

17. We may apply these formulae to develope any positive 
integral power of /t in a series of zonal harmonics, as we 
proceed to shew. 

Suppose that m is a positive integer, and that /a*" is de- 
veloped in such a series, the coeflScient of i^ being C^ so 
that 

then, multiplying both sides of this equation by P, and inte- 
grating between the limits — 1 and 1, all the terms on the 

right-hand side will disappear except I C^I^^dfl, which will 

2 ' 

become equal to q. .. (7,. 

2/ -4-1 r* 
Hence ^* = ^^^J /*"-P*^/*' 

which is equal to 0, if m + 1 be odd. Hence no terms appear 
unless m + « be even. In this case we have 



«(2t + l)[ pTPtdfi. 

J Ik 



+ 7 , -- ... P.+ ^I>, 
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Hence the formula just investigated gives 
^ _ /o . . 1 N (m-t + 2)(m-i+4)...(m-l) 
^*"^'^*'*"^(m+i + l)(m + «-l)...(m + 4)(m + 2) 

if i be odd, and 

r - W J- -W (m-^ + 2)(m-/+4)...m 

o, - ^^ii- 1; {m + i+ 1) (m+ »-l) ... (m + 3) (m+ 1) 

if t be even. 

Therefore if m be odd, 

m^/a i\ 2.4.6 ... (m — 1) p 

/i, -C^^ + l;(2m+l)(2m-l)...(w + 4)(»i + 2)^«"^- 

m— 1 

(m + 4)(m + 2)"« ' m 

If m be even, 

/* = (^ + i;^2m+ 1) (2m-l) ... (m+3) (m + l)'^""*' "* 

I 5 "^ p , 1 p 

^ (m + 3)(m + l) »^m+l •" 

Hence, putting for m successively 0, 1, 2 ... 10, we get 
/ = -?«. 

„ t _ K _2_ pii:p_2plp 

it'' = 7=^P+-P 

'^ 7.5 » 5 » 

=?p+?p 

5 » 5 >' 

35 * 7 «'*'5 •* 
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'^ 11.9.7 "^9.7 »^7' 
63 '^9 »^7 " 

231 « 77 * 21 * 7 *^ 

i.» — 1 g 2.4.6 p , T 1 ^-^ p X 7 ^ P J. ^ P 
'^ ~^^15.13.11.9"'''*'^^13.11.9''»^'ir:9^»'^9 " 

. _iipj.Ap j.1*P4.1p 
~429 ' 39 * 33 » 3 " 

»_ -ihr ^.4.6.8 p ,q 4.6.8 p 
'*~;^' 17. 15. 13. 11.9-^' "^""^ 15. 13. 11. 9"^' 

j.q ^-^ P 4.K_§_ p J.1 p 
^13.11.9^^ 11.9 «^9'" 

_128^ 64 48 40 1 

~ 6435 » y 495 ' ^ 143 "^^ ^ 99 » ^ 9 »' 

„<•_ iQ 2.4.6.8 p , ,5 '*-6-8 o 

'^ ~ 19.17.15.13.11 •^17.15.13.11^ 

a. n _Al?__ p J. 17 _§_ P 4. q J. p 
■*■ 15.13.11 •■^'13.11 »■*' 11 " 

128 p^i92p 16 ^ J^ 



12155 "^2431^7^65 '^143 '^ll 

2.4.6.8 .10 p -y 4.6.8.10 p 

'^ ~ 21 . 19 . 17 . 15 . 13 . 11 ""^ 19.17.15.13.11 ' 

,.,, 6.8.10 p . ^ 8.10 p 16 p.lp 

''""l7.15.13.11"^» "^^15.13.11 ■^*"*'^l3ai^«"^ 11 "^0 

_ 256 p 128 32 48 p 50^ Ip 

~ 46189"»"^2717 ' 187 V"*" 143 '"^ 143 ^"^ 11 •* 
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18. Any zonal harmonic P, may be expressed in a finite 
series of cosines of multiples of 0, these multiples being 
i0,{i-2)0.... Thus 

therefore, writing cos for /a, and observing that 

1 -2 cos^A + A»= (1 - Ae^^^O (1 - Ae-^^«), 
we obtain 

or 



(H.|A€>^-1« + |^A*€V-12^ + ... 



-^^^S^*"^--"--) 



(l + 1 Ae->^» + 1^^ AVV^w + ... 



-^^?s^ »•«—-••■) 



whence, equating coefficients of h\ 

■P«= a . o- ' 2cosid + ir-A )ir= — 5ro2cos(»-2) ^ 

* 2. 4.. .2* 2.4... {2i-2)2 ^ ' 

^1.3...(2t-5)1.3_ .. ,.^^ 

+ 2747:7(271:4) 2:42<=°«(»-*)^+-- 

the last term being •! ' *' — ; — i- if i be even, and 

I ^ • tB , • t % J 

. 2.4...(i + l)2.4...(i-l)^ ^°^^'^^^""'^^- 
19, Let us next proceed to investigate the value of 

I P^ cos m0 sin d d0* 
Jo * 
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This might be done, by direct integration, from the above 
expression. Or we may proceed as follows. 

The above value of P^ when multiplied by cos mO sin 

(that is by ^ {sin (m+ 1) ^ — sin (m — 1) ^}) will consist of a 

series of sines of angles of the form {i — 2ri + (m + 1)} 0, that 
is of even or odd multiples of 0, as i + m is odd or even. 
Therefore, when integrated between the limits and ir it 
will vanish, if i + m be odd. We may therefore limit our- 
selves to the case in which i + m is even. 

Again, since cos m0 can be expressed in a series of powers 
of cos 0, and the highest power involved in such an expression 
is cos"*^, it follows that the highest zonal harmonic in the 

development of cos m0 will be P^. Hence / P^ cos m0 sin d0 

Jo 

will be = 0, if m be less than i. 

Now, writing 

P, = C; cos ;5+ (7^ cos (t- 2) ^ + ... 

we see that P, cos m0 sin d0 will consist of a series of sines 
of angles of the forms (m + i+V) 0, {m + i — l) ... down to 
(m — t— 1)^, there being no term involving m0, since the 
coefficient of such a term must be zero. Hence 



Jo 



P^ cos m0 sin d0j 



will consist of a series of fractions whose denominators in- 
volve the factors m + i+l,m + i^l ... m-^i—l respectively. 
Therefore when reduced to a common denominator, the result 
will involve in its denominator the factor 

(m + i+ 1) {m + i—l) ... (m + 1) (m — 1) ... (wi — i — 1) 
if m be even, and 

(m + e + 1) (w + t-1) .-. (m + 2) (m-2) ... (w-t- 1) 
if m be odd. 

For the numerator we may observe that since 

cos m sin 0d0 



S?' 
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Vanishes if m be less than t, it must involve the factors 
in — (i — 2), m— (i — 4) ... 7n + (i — 2), and that it does not 
change sign with m. Hence it will involve the factor 

{m-(t-.2)}{m-(i-4)}..,(w-2)m«(w + 2)...(w + f-2) 

if m be even, and 

{w-(;-2)}{m-.(i-4)}... (w-.l)(m + l)..-(m + ;-2) 

if m be odd. 

To determine the factor independent of m, we may pro- 
ceed as follows : 

P, = (7, cos 10 + G^ cos (i- 2) 5 + ... ; 
/. P^cos w^= X (7, {cos (w + i) ^ + cos (m — t) ^} 

+ 2 ^*-« f^^^ (^ ■*■ *" ^^ ^ + ^^^ (w - i + 2) ^} + ... ; 

1 
.'. P,cosm^sin^ = J C^fsin (w + i + 1) ^ — sin (??» + i— 1) 5 

H- sin (m — « +1) ^ — sin (m — i — 1) ^} 
+ T C't-a {sin (w + * — 1) ^ — sin (m + 1 - 3) 



•' 



+ sin {m — / + 3) 5 — sin (wi — t + 1) ^} + . . . ; 
cos mO sin ^ c?5 



= A|_i l_+— 1 ?_l 

2 \m + i-l 7/i4-i— 3 m-t + 3 m — i+lj 
~ ^' 1 m' - (i + 1/ ^ m« - (t - 1)*; 
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Now, when m is very large as compared with %, this be^ 
comes 



-2 



C«+ Cf^+ ... 



V 2 y 



since G^ + (7^, + ... = 1, as may be seen by putting ^ = 0. 

Hence I P. cos mO sin ^ d5 tends to the limit . ^ as m 

./o / mr' 

is indefinitely increased. 

The value of the factor involving m has been shewn 
above to be 

{m -- {% " 2)} {m - (t- 4)} ... (m~ 2) m'(m + 2) ... (;n 4-t - 2) 
{m - (1 + 1)1 {^ - (* - i)l ••• (w - 1) (w + 1) ... (m + * + 1) 
if m be even, and 

{^,(^'.2)}{m-(t-4)}...(m-l)(m + l)... (m + t-2) 
{m--(i+l)}{m-.(i-l)}... (w--2)(w + 2)... (w + i+1) 

if m be odd. 

Each of these factors contains in its numerator two factors 
less than in its denominator. It approaches, therefore, when 

m is indefinitely increased, to the value -^ . Hence 



/, 



P^ cos mO sin 6 dd 





{m--(t-2)]{m~(t-4)]...(m~2)m»(m+2)...{m+(t~2)} 

" {wi-(*+l)}{m-(i--l)}...(w-l)(w + l)...lm+(i+l)} 

if m and i be even, and 

._ __ { m-(t-2)}{m-(t-4)]...(m-l)(m+l)...{m+(t-2)} 

"" {w-.(i + l)}{m-(t-l)}...(m-2)(m+2)...{m+(i + l)] 

if m and i be odd* 

In each of these expressions i may be any integer such 
that m — t is even, i being not greater than m. Hence they 
will always be negative, except when % is equal to m. 
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20. We may apply these expressions to develop cos mO 
in a series of zonal harmonics. > 

Assume 

Multiply by Pi sin 6, and integrate between the limits 
and IT, and we get 

^ {m- (t - 2)} {m - (^• - 4)} ... {m + {{- 2)} 2_ 

{»i- (» + 1)] {to - («■- 1)1 ... [m + (t + 1)} ii+ 1 '• 

Hence 

^^' "*■ ^^ {m - (;+ 1)1 {m - (2 - 1)1 ... {w + (i + 1)1 

Hence, putting w successively = 0, 1, 2, ... 10, 
cosO^ = P,; 
cos 6 = P^; 

3 > 3 ^' 

6 » 6 *' 



3.5 • 



« 64 p _ 16 p _ J_ p . 
35 * 21 « 15 *' 



F.H. 
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KO— -n 2.4.6.8 p fj 4.6 p 



~ 63 » 9 » 7 " 
.n«fi/)- iq 2. 4. 6'. 8. 10 p 

4.6'.8 6' p 1 

1.3.5.7.9.11 * 3.5.7.9 * 5.7 ' 

= ^P_38ip_lp_lp. 
231 ' 385 * 21 * 35 •' 

^a ir 2.4.6.8.10.12 „ 

cos 7^ = - 16 -j-j-g-g-g-Jj-Jg-^ P, 

_ 4.6.8.10 p y 6.8 p 3^ p 

1.3.5.9.11.13 » '3.5.9.11 » 5.9 ' 

1024 128 112 J^ 

■'•7 TT7^6 AOK-'^J 1K'*1' 



cos85 = -17 



429 ' 117 " 495 » 15 
2.4.6.8M0.12.14 



-1.1.3.5.7.9.11.13.15.17 ' 



4.6.8M0.12 p 6.8M0 p 

1.3.5.7.9.11.13.15 ' 3.5.7.9.11.13 * 

K ^' P L_P 

6.7.9.11 » 7.9 • 

16384 p 4096 p 256 p 64 p 1 p 
~ 6436 ' 3465 • 1001 * 693 » 63 • ' 
coaQg- 19 2--t-6-8.10.12.14.16 ' 

4.6.8.10.12.14 p 6.8.10.12 

1.3.6.7.11.13.15.17 ' 3.5.7.11.13.15 » 

8.10 p 3 1 p 
'6.7.11.13 » 7.11 » 
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_ 32768 3072 „ 128 „ 16 „ 3 „ 
"12155 • 2431 ' 455 • 143 » 77 " 

._ 2.4.6.8.10M2.14.16.18 „ 

cos lot/ ^i_i. 1.3. 5. 7. 9. 11. 13. 15. 17. 19. 21^" 

y 4. 6. 8. lO*. 12. 14. 16 p 
1.3.5.7.9.11.13.15.17.19 » 

6.8. 10*. 12.14 p 8.10M2 p 

3.5.7.9.11.13.15.17 • 5.7.9.11.13.15^ 

5 10' p Lp 

. 7.9.11.13 » 9.11 • 

_ 131072 _ 32768 p _ _512 ^28 » _ ^^ p 

~ 46189 " 24453 ' 1683 • 1001 * 9009 » 

99 •■ 

21. The present will be a convenient opportunity for 
investigating the development of sin in a series of zonal 

harmonics. Since sin ^ = (1 — /*')i it will be seen that the 
series must be infinite, and that no zonal harmonic of an odd 
order can enter. Assume then 

sind= C,P,+ C^P^ + ... + CiPi+ ... 

t being any even integer. 

Multiplying by P«, and integrating with respect to /* 
between the limits — 1 and + 1, we get 



/ 



1 2 



Jo 



2 .,0 

supposing P{ expressed in terms of the cosines of and its 
n^ultiplea 

^?i+irp^(l.CQs2d)dft 

or— —2* 
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Hence, putting i = 0, 

13 

Putting t = % and observing that JP, = j + t cos 2^> 

„ 5 f' (1 + 3 cos 20) (1 - cos W) J. 
^*"4Jo 4 '^ 

= ^ r jl + 2 cos 2^ - 1 (1 + cos 4^)| <f« 

5 
= "32''- 

For values of i exceeding 2, we observe, that if we write 
for Pi the expression investigated in Art. 18, the only part 

of the expression | P* (1 — cos 20) dO which does not vanish 

Jo 

will arise either from the terms in P< which involve cos 25, or 
from those which are independent of 6, We have therefore 

C, = 2^41 p 1-1.3^4+ 1) |4-f:| 2 cos 2^ 
4 Jo l2.4... (i + 2) 2.4 ... (i — 2) 

_ 2e + l 1.3...(t-l)1.3...(i-3) 
4 '2.4... i 2.4... (t-2) 

{' (irl + *±_| 2 COS 2^) (1 - cos 20) d0 

_ 2t + 1 1 .3 ... (t- 1) 1 .3 ... (t- 3 ) /t-1 t + l\ 
~ 4 2.4... i 2A...{i-2)'^\ i t + 2/ 
2t + 11.3... {i-V) 1.3... (i-3) 



= — w 



2 2.4...«(t + 2)2.4...(t-2)* 



Hence Bin^ = T Po — qo ■'*«"••• 

(2t + l)7rl.3... (t-1) 1.3... (ir-3) p 
2 2.4...t(i + 2)2.4... (»-2)t-^*~"* 
i being any even integer. 
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dP 

22. It will be seen that -~ , being a rational and integral 

function of /i»*"\ /i*'*'..., must be expressible in terms of 
Pi-xi -^i-s''' T^ determine this expression, assume 

dP 

then multiplying by P,», and integrating with respect to //. 
from — 1 to + 1, 

And /p„g'rf;* = P.P,-/p.^ci/*. 

Now, since i>m, 

since either i» or % mnst be odd, and therefore either P„ or 
P, = — 1, when /t = — 1 ; 



2^ 
2to + 1 



.-. 2 = ^ . , C^ or C,= 2»i + 1; 



/. ^* = (2i - 1) P^, + (2i - 6) P... + (2» - 9) P^, + . . . 

23. From this equation we deduce 

P.~P^— (2»*-l)fP^,dl*. 

the limits ii and 1 being taken, in order that P, — P^.g may 
be equal to at the superior limit. 
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Now, recurring to the fundamental equation for a zonal 
harmonic, we see that 

• P—P = — ^*~^ (\ _ u*\ ^^*^ 

-(5n'^)a-'")f'- 

24. We have already seen that / P,P«, d/A == 0, i and m 
being different positive integers. Suppose now that it is 
required to find the value of | P^IP^ ^/^» 

We have already seen (Art. 10) that 

iu. (t-m) (i + m+1) 



And, from above, 



^^ *^'^ — 2m +1 ^^"« ^-'^* 

» (* + 1) -p (-p p n1 
~ 2» + l " ^ *** ■^'-i-'J- 



ZONAL HARMONICS. 39 

25. We will next proceed to give two modes of ex- 
pressing Zonal Harmonics, by means of Definite Integrals. 
The two expressions are as follows : 






dS^ 



{/^±0iA»-l)4cos&}'*^' 

IT Jo 

These we proceed to establish. 
Consider the equation 



-r- 



d^ 



ftcosS^ (a'-6")4' 

The only limitation upon the quantities denoted by a 
and b in this equation is that V should not be greater than 
a\ For, if 6' be not greater than a\ cos & cannot become 

equal to j- while Sr increases from to tt, and therefore the 

expression under the integral sign cannot become infinite. 

Supposing then that we write z for a, and V— 1 p for 6, 
we get 

1 r» d^ ^ 1 

irJo « — V^pcosSr («'+f)')*' 

We may remark, in passing, that 



Jo « — V— IpCOsSr Jo ^ + V— IpCOSSr 



-r 



zddr 



«" + f)' cos' St * 

and is therefore wholly real. 

Supposing that p*=^a? + j^, and that iaj'+y* + «* = r*, we 
thus obtain 

If' d^ ^ 1 

^Jo ^-V^pCOSSr ^* 
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Diflferentiate t times with respect to z^ and there results 

i^f' a^ _d'i p 

Hence P=— J- ^^i_. J!. f ^ 
^1*^ f' d^ 

In this, write fir for z, and (1 — fj?)^ r for p, and we get 
1 r» ^ 

^*"^Jo {M-(/i''-i)*cos^r'' 

which, writing tt— Sr for Sr, gives 



1 /•» d 

"^* " ^ Jo Im^O^ 1)* cos^l'^^ • 

26. Again, we have 

1 1 f^ dyjr 

(a*— 6')* IT Jo a — bcosy^' 

In this write 1 —fihtoT a, and + (jj? — !)^ hior 6, which is 

admissible for all values of h from up to /a — (/*' — 1)* and 
we obtain, since a' — 6* becomes 1 — 2/aA + h\ 

1 1 r* dy^ 

(1-2M + A')*"^Jo1-M+0m;*-1)* 



Acos*^ 
dyfr 

TT. 



If' djr . 

TtJo l-{/A±(/Lt*— l)*COS'^|r} a' 



= ^r^^[l + {/*±(/^'-l)*COS^}A+... 

+ {a* ± 0*' - 1)* cos fyh* +...]• 
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Hence, equating coefficients of A', 

The equality of the two expressions thus obtained for P^ is 
in harmony with the fact to which attention has already 
been directed, that the value of P^ is unaltered if — (i + 1) 
be written for i. 

27. The equality of the two definite integrals which 
thus present themselves may be illustrated by the following 
geometrical considerations. 




Let be the centre of a circle, radius a, C any point 
within the circle, PCQ any chord drawn through C, and let 
0(7=6, GOP^%COQ-^'>^. Then(7P»=a' + i»-2aJcos&, 
C^ = a' + 6' — 2a6 cos -^^ Hence 

(a«+6*-2a6cos^)(a* + y-2aicos^) = (a»-6y; 

sin^c?^ sin'^rf^ _^ 

Va« + 6»-2a6cos^"*'a'+i'-2a6co8^" • 
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Again, since the angles OP (7, OQO sue equal to one 
another, 

sin _ sin OPC _ sin OQG _ sin yfr 
CP OU " 00 " CQ' 

. sinS- sin-^^ 



"(a=+6»-2a6cos&)*""(a* + 6'-2a6cos^)4' 

whence . • ^^^ r = 0. 

(a*+6'-.2a6cos^)4^(a' + 6'-2a6cosf)4 

^^ / 2^x, o ^ ST;r4=(a' + i"-2a5costrJ; 

(a* + 6*— 2aocosa^)^^ ^ 

.-. (a« - &^)«^^ , . ,, f\ ^T^= - (a' + &•- 2a6cos f)' c?^. 

"^ ' (a'*+6' — 2aJcos&)*** ^ r/ r 

In this, write a' + 6' = /i», 2db =?(/*'— 1)^, which gives 
a^'-b^^ 1, and we get 

We also see, by reference to the figure, that as Sr in- 
creases from to w, -^ diminishes from tt to 0. Hence 

J {/A + (/Lt* — 1)^ COS &}**^ Jo 

28. From the last definite integral, we may obtain an ex- 
pansion of P| in terms of cos 9 and sin 9. Putting fM = cos 9, 
we get 

+ {cos 9 — V^ cos y^ sin ^}*] d^ 
= - r{(cos^y-%=^cos«t(cos^*-'(sin^)«+... 

TT J Q ± • Z 
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+ (_!)« »(t'-l)-.(t-2^ + l) (pj,g ^« (^^ 0y^ (gjjj ^^ 

X • ^ • • • ^772) 

1VT f'/ .\«mj. (2m-l)(2OT-3)...l 

Now J^(co8frJV-=^ ^ 2m(2m-2)...2 ' 

A A »(»'-l)-(*-2"t + l) (2m-l)(2TO-3)...l 

1.2.. .2m 2TO(2m-2)...2 

_ ^(^•-l)...(^•-2m+l) . 
(2.4...2m)' ' 

.-. P, = (cos 0)' - :^^^ (cos ^*^ (sin ^* + . . . 



CHAPTER III. 

APPLICATION OF 20NAL HARMONICS TO THE THEORY OP 
ATTRACTION. REPRESENTATION OF DISCONTINUOUS 

FUNCTIONS BY SERIES OF ZONAL HARMONICS. 

1. We shall, in this chapter, give some applications of 
Zonal Harmonics to the determination of the potential of a 
solid of revolution, symmetrical about an axis. When the 
value of this potential, at every point of the axis, is known, 
we can obtain, by means of these functions, an expression 
for the potential at any point which can be reached from 
the axk Without pa.sing\Lugh the attra^ing ma«a. 

The simplest case of this kind is that in which the 
attracting mass is an uniform circular wire, of indefinitely 
small transverse section. 

Let c be the radius of such a wire, p its density, k its 
transverse section. Then its mass, M, will be equal to 2'7rpck, 
and if its centre be taken as the origin, its potential at any 

M 

point of its axis, distant z from its centre, will be — r« 

Now, this expression may be developed into either of the 
following series : 



M 
z 



f 2?"*"2.4;5* '"^^ ^ 2.4...2i ^ai + -|-v^;. 

We must employ the series (1) if « be less than c, or 
if the attracted point lie within the sphere of which tlie ring 
is a great circle, and the series (2) if « be greater than c, 
or if the attracted point lie without this sphere. 
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Now, take any point whose distance from the centre is r, 
and let the inclination of this distance to the axis of the 
ring be B. In accordance with the notation already em- 
ployed, let cos ^ 3B /A. Then, the potential at this point will 
be given by one of the following series : 



1.3.5...(2t-1) r" 1 



2.4.G...2i 






r ( 



7^ • 2.4 



^^ ^^ 2.4.6...2i ^«'^ + --|-^^^ 

For each of these expressions, when substituted for F, 
satisfies the equation V*F=0, and they become respectively 
* equal to (1) and (2) when Q is put = 0, and consequently 
?• = «. The expression (2') also vanishes when r is infinitely 
great, and must therefore be employed for values of r greater 
than c, while (!') becomes equal to (2') when r = Cy and will 
therefore denote the required potential for all values of r 
less than c. 

These expressions may be reduced to other forms by 
means of the expressions investigated in Chap. 2, Art. 26, viz. 

P, = i r(/i^ + ^;7^cos^yd&, 
or P, = i rO* + V;?=Tcosfr'^^^cZt. 

TT Jo 

Substitute the first of these in (!') and (observing that 
yi/r = z) we see that it assumes the form 



-M;r'f i {^+(g'-0*cos^}' 

TTcJot 2 ^ 



. 1.3 {^-Kig' -O^cos^r Lycv 
^ 274 ? -j '^^' 
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whiclL is equivalent to 

Mr d^ 

The substitution of the last form of P^ in the series (2') 
brings it into the fonn 

Mr ( 1 1 c' 

which is equivalent to 

Mr ^ 

2. Suppose next that the attracting mass is a hollow shell 
of uniform density, whose exterior and interior bounding 
surfaces are both surfaces of revolution, their common axis 
being the axis of z. Let the origin be taken within the 
interior bounding surface ; and suppose the potential, at any 
point of the axis within this surface, to be 

AQ + A^z + Aji^+ ... +A/+ ... 

Then the potential at any point lying within the inner 
bounding surface will be 

AJ',+A,P,r + A^y+... + A,P/ + ... 

For this expression, when substituted for F, satisfies the 
equation V^F=0; it also agrees with the given value of 
the potential for every point of the axis, lying within the 
inner bounding surface, and does not become infinite at any 
point within that surface. 

Again, suppose the potential at any point of the axis 
without the outer bounding surface to be 

JDq xf. ^o ^i 
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Then the potential at any point lying without the outer 
bounding surface will be 

^0^0. Ml. M. .^^^^4. 

For this expression, when substituted for V, satisfies the 
equation V* F= ; it also agrees with the given value of the 
potential for every point of the axis, lying without the outer 
bounding surface, and it does not become infinite at any 
point within that surface. 

By the introduction of the expressions for zonal har- 
monics in the form of definite integrals, it will be found that 
if the value of either of these potentials for any point in the 
axis be denoted by ^ (z), the corresponding value for any 
other point, which can be reached without passing through 
any portion of the attracting mass, will be 



i r^{i5 + («*-04cos^}d^. 



3. We may next shew how to obtain, in terms of a series of 
zonal harmonics, an expression for the solid angle subtended 
by a circle at any point. We must first prove the following 
theorem. 

The solid angle, svhtended hy a closed plane curve at any 
point, is proportional to the component attraQtion perpendicular 
to the plane of the curve, exercised upon the point by a lamina, 
of uniform density and thickness, bounded by the closed plane 
curve. 

For, if dS be any element of such a lamina, r its distance 
from the attracted point, the inclination of r to the line 
perpendicular to the plane of the lamina, the elementary 
solid angle subtended by d8 at the point will be 

dS cos 

• r- ; 

And the component attraction of the element of the 

lamina corresponding to dS in the direction perpendicular 

to its plane wiU be 

7 dS ^ 
pk —^ cos u, 



48 APPLICATION OF ZONAL HAEMONICS 

p being the density of the lamina, Jc its thickness. Hence, 
for this element, the component attraction is to the solid 
angle as pk to 1, and the same relation holding for every 
element of the lamina, we see that the component attraction 
of the whole lamina is to the solid angle subtended by the 
whole curve as pk to 1. 

Now, if the plane of {cy be taken parallel to the plane 
of the lamina, and V be the potential of the lamina, its 
component attraction perpendicular to its plane will be 

— -7—. Now since Fis a potential we have V'F=0, whence 

-7-V*F=0, or ^*{";7~)=0. Hence --r- is itself a potential, 

and satisfies all the analytical conditions to which a potential 
is subject. It follows that, if the solid angle subtended by 
a closed plane curve at any point (x, y, z) be denoted by 
CO, o) will be a function of x, y, z, satisfying the equation 
V^o) = 0. Hence, if the closed plane curve be a circle it 
follows that the magnitude of the solid angle which it sub- 
tends at any point may be obtained by first determining 
the solid angle which it subtends at any point of a line 
drawn through its centre perpendicular to its plane, and 
then deducing the general expression by the employment 
of zonal harmonics. 

Now let be the centre of the circle, Q any point on the 
line drawn through perpendicular to the plane of the 
circle, E any point in the circumference of the circle. With 
centre Q, and radius QO, describe a circle, cuttiDg QE in L. 
From L draw LN, perpendicular to Q 0. 

Let 0^=c, OQ=-z. 

Then^i = (c«+^)4.^, ON^ -^^i {{(? + z^)i ^ z] 



(c' + z')i' 
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And the solid angle subtended by the circle at 

A ON 



= 27rll '—\. 

1 (c-+;^)ij 



To obtain the general expression for the solid angle sub- 
tended at any point, distant r from the centre, we first 
develope this expression in a converging series, proceeding 
by powers of z. This will be 

^X c'^2c» 2.4c' ^ ' 2A...2i ^' + -| 

if ^ be less than c. and 

9^ il ^^ _ 111 £* J. f iv l-3-.(2t-l) c'', ' 

if 3 be greater than c. 

Hence, by similar reasoning to that already employed, 
we get, for the solid angle subtended at a point distant r 
from the centre, 

^'^r* +2 c' 2.4 c" ^"' 

^ '' 2.4...2i c"*' ■*"*••) 
if r be less than c, and 

[2 r" 2.4 r* ■*""• ^ ^ 2.4...2i r''* "*"••; 
if r be greater than c. 

4. We may deduce from this, expressions for the potential 
of a circular lamina, of uniform thickness and density, at an 
external point. For we see that, if V be the potential of 
such a lamina, Ic its thickness, and p its density, we have for 
a point on the axis, 



-S=''*-2-{i-(^i}' 



F. H. 
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whence F= lirph {(c* + z^)^ — z] 

if M be the mass of the lamina. 

Expanding this in a converging series, we get 






. 1«» 1.1 2* 1.1.3/ 



2c 2. 4 0" ^2. 4. 6c' '"' 

1.1.3...(2t-3) g'' 1 

"^"■^^ 2.4.6...2t c'^«"*'-j 

ff 

if i2? be less than c, and 

^"c^X^z iAz'^'^^A.Qz' •"• 

1.1.3...(2i~3) c'' 1 

^^"•^^ 2.4. 6.. .2i 2«^i"^-j 

if z be greater than c. 

Hence we obtain the following expressions for the po- 
tential of an uniform circular lamina at a point distant r 
from the centre of the lamina : 



r=f{p,c-p,r+l 



1 py 1.1 py 



a ~ • • • 



c 2.4 c" 

. ■,., 1.1.3...(2t-3)P^r" , \ 
^ ^' 2.4.6...2i c"-* ■^••■•j 

if r be less than c, and 

.r_M\\ P/ 1.1 P,c« 1 . 1 .3 P/ 
*^~c't2 r 2.4 r» "^2.4.6 i* '" 

, ,,. 1.1.3...(2t-3)P^c" . \ 
~^~^'' 2.4.6...2t r«-^ ■^••••J 

if r be greater than c. 
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It may be shewn that the solid angle may be expressed 
in the form 

2^-2 r ^ + (^'-0\eos^ 

^ [c'+ {^ + is? - r')i cos ey]^ 

and the potential of the lamina in the fonn 



e 



^r[c'+{^ + (;s»-0*cos^nic?^-^. 



5. As another example, let it be required to determine 
the potential of a solid sphere, whose density varies inversely 
as the fifth power of the distance from a given external point 
at any point of its mass. 

It is proved by the method of inversion (see Thomson 

and Tait's Natural Philosophy, Vol. 1, Art. 518) that the 

M 
potential at any external point P' will be equal to- p, , 0' 

being the image of in the surface of the sphere, and M 
the mass of the sphere. We shall avail ourselves of this 
result to determine the potential at a given internal point. 

Let G be the centre of the sphere, the given external 
point. Join CO, and let it cut the surface of the sphere in A, 
and in CA take a point (X, such that CO.CO^ CA\ Then 
(y is the image of 0. 

Let P be any point in the body of the sphere, then we 
wish to find the potential of the sphere at P. 

Take as pole, and OG as prime radius, let OP = r, 
POG^e. Alsolet(7^ = a, (70 = c. 

Let the density of the sphere at its centre be p, then its 
density at P will be p -5 . Hence 



Jf = 27r JJp p 7^ sin ^ dr d^, 



4-^2 
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the limits of r being the two values of r which satisfy the 
equation of the surface of the sphere, viz. 

r* + c" — 2cr cos ^ = a*, 
and those of 6 being and sin"* - . 

Hence, if r^, r^ be the two limiting values of r, we have 



^ 1 1 2c COS tf /I 1\ 

Now -^ ^=-g ^ 1. 

., 1 1 2c cos ^ 

Also r + r = -:» % » 



1 



H 8 9 



. /I l\^g {c'co8«g-(c*-a*))4 
\r, Tj/ c* — a* 



_- fa« -c' sin' g)^ 



.-. i/ = — i^ ij a . -^ 2 COS sm ^ (a'- c' sin'^)^ 

= ^^^, J'*"''^ cos ^ sin ^ (o' - c» sin* ^ » rf^ 

4 TTjOC* 



(^5 






Now, if F be the potential at P, we have (see Chap. i. 
Art. 1) 

^{rV) ,1 d f . adV\ 4'7rpc' ' 



1' — j„ — + r-- 






dr' • sine ddV deJ r* 
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r 

This is satisfied by F= — ^ — =-. 

Assume then, as the complete solution of the equation, 



(j/+^)p.+.... 



It remains to determine the coefficients A^, A^.,.Ai..,B^, 

B^...B^i so that this expression may not become infinite for 

any value of r corresponding to a point within the sphere, 

and that at any point P on the surface of the sphere it may 

M 
be equal to -jf-py where CfP : OP :: a:c, and therefore, at 

the surface, 

tt_Mg 1 _4 TTpc^a* 

And, at the surface, we have 
r*- 2crcos^ + c'^ — a' = 0; 

1 _ 1 /I 2ccosg \ 

cr — a \ r rj 

+ (^ir+ -^'1 P, + ... identically. 
,T 4 -TrpcV /'2 -TTpC* , „ \ D 

and -Bj, -B„...B,... J.^, .4j...-4, all = 0. 
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Hence since P, = 1, 



2 nrp^ ( 2a* ,) 
•~3c*-o'l(c'-aO I 

-^ ^P^ r^/i' r'\ 



whence we obtain, as the expression for the potential at any 
internal point, 

6. We shall next proceed to establish the proposition that 
if the density of a spherical shell, of indefinitely small thick- 
ness, he a zonal surface harmonic, its potential at any internal 
point will he proportional to the corresponding solid har- 
monic of positive degree, and its potential at any eocternal 
point mil he proportional to the corresponding solid harmonic 
of negative degree. 

Take the centre of the sphere as origin, and the axis of 
the system of zonal harmonics as the axis of z. Let h be the 
radius of the sphere, hh its thickness, U its volume, so that 
U= iirb^Bh. Let GP^ be the density of the sphere, P, being 
the zonal surface harmonic of the degree i, and C any con- 
stant. 

Draw two planes cutting the sphere perpendicular to the 
axis of z, at distances from the centre equal to f, f + df 
respectively. The volume of the strip of the sphere inter- 

cepted between these planes will he ~ U, and its mass will be 

CP,§U. 

Now ^==hfif hence d^ = hdfi, and this mass becomes 
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Hence the potential of this strip at a point on the axis of z^ 
distant z from the centre, will be 

GU P| J 

2 (3« + 6» - 2bz(i)^ ^' 



which may be expanded into 

CU 
2b 



P. fP,+ P^ 1+ ... + P, ^+ ....) d(li£3<b, 

and ^^(p^ + F^^^+... + p^^^ + ...^d^i£z>b. 

To obtain the potential of the whole shell, we must inte- 
grate these expressions with respect to fi between the limits 
— 1 and + 1. HeDce by the fundamental property of Zonal 
Harmonic^, proved in Chap. II, Art. 10, we get for the po- 
tential of the whole shell 

^. -. jTfi at an internal point, 

9 'jIT "<+i ^^ ^^ external point. 

From these expressions for the potential at a point on 
the axis we deduce, by the method of Art. 1 of the present 
Chapter, the following expressions for the potential at any 
point whatever : 

CU Pr* 
r, = ,-r^ % at an internal point, 

From hence we deduce the following expressions for the 
normal component of the attraction on the point. 

Normal component of the attraction on an internal point, 
measured towards the centre of the sphere, 

~~ dr ~~2i+l^'^ 6'«~* 



i 
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Normal component of the attraction on an external point, 
measured towards the sphere, 

In the immediate neighbourhood of the sphere, where r is 
indefinitely nearly equal to J, these normal component at- 
tractions become respectivdy 

and their diflFerence is therefore 

p 

And writing for TJ its value, 47rJ*S5, this expression be- 
comes 

47rS6 . CP,. 

Or, the density may be obtained by dividing the alge- 
braic sum of the nonnal component attractions on two points, 
one external and the other internal, indefinitely near the 
sphere, and situated on the same normal, by 47r x thickness 
of the shell 

7. It follows from this that if the density of a spherical 
shell be expressed by the series 

CiP,+ C7,P,+ C7^,+ ... + 0,P, + ..., 

G^y C^y (7, ... <7, ... being any constants, its potential (V^) at 
an internal point will be 

\ b ^3 b* ^5~P~^'"^2t + l V*' ^""J 
and its potential ( F^ at an external point will be 

In the last two Articles, by the word ''density** is meant * 
"volume density," i.e. the mass of an indefinitely small 
element of the attracting sphere, divided by the volume of 



iirV 
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the same element. The product of the volume density of 
any element of the shell, into the thickness of the shell in 
the neighbourhood of that element, is called "surface den- 
sity." We see from the above that, if the surface density 
be expressed by the series 

the potentials at an internal and an external point will seve- 
rally be 

,/^ 1^^ 1^^ 1 cr.fiy N 

This variation in surface density may be obtained either 
by combining a variable volume density with an uniform 
thickness, as we have supposed, or by combining a variable 
thickness with a uniform volume density, or by varying both 
thickness and density. 

8. We have seen, in Chap, ii., that any positive integral 
power of iJb, and therefore of course any rational integral 
function of fi, may be expressed by a finite series of zonal 
harmonics. It follows, therefore, that we can determine the 
potential of a spherical shell, whose density is any rational 
integral function of /a. 

Suppose, for instance, we have a shell whose density 
varies as the square of the distance from a diametral plane. 
Taking this plane as that of xy^ the density may be ex- 

pressed by p^y oi pjt. We have seen (Chap. n. Art. 20) 

that 

/^» = |(1+2PJ. 

Hence, by the result of the last Article, the potential 
will be 

p -^ (r + K -hr) at an internal point. 
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p o^ ( - + e — V" J 8»t an external point ; 
or, since Fjr = -^--^ — r = — ^ — , we obtain 
p ^ f r + ^ — o — j for the potential at an internal point, 

p -o 1~ + t("";3 ■* — »")r ^"^^ *^^* ^* ^^ external point. 

9. As an example of the case in which the density is re- 
presented by an infinite series of zonal harmonics, suppose we 
wish to investigate the potential of a spherical shell, whose 
density varies as the distance from a diameter. Taking this 
diameter as the axis of z, the density will be represented 

by p sin d, or p (1 —jjf)^. We have investigated in Chap. II. 
Art. 21, the expansion of sin d in an infinite series of zonal 
harmonics. Employing this expansion, we shall obtain for 
the potential 



TTpU 

2T 



2 • 16^6^ ••• 2.4..i(i+2)-2.4...(i-2)* 'J* '"y 



or 



^ ttH^^Lp^I^ i»3...(t-i) i.3...(i-3) p y ) 

2^^|2r 16 V ••• 2.4...t(i+2)'2A..(i-2)i V^^ "']; 

according as the attracted point is internal or external to the 
spherical shell, t being any even integer. All these expres- 
sions may be obtained in terms of surface density, by writing, 
instead o{ pU, 47rcV. 

10. We may next proceed to shew how the potential of 
a spherical shell of finite thickness, whose density is any solid 
zonal harmonic, may be determined. Suppose, for instance, 
that we have a shell of external radius a, and internal radius 
a', whose density, at the distance c from the centre, is 

■ri PiC*, h being any line of constant length. 

Dividing the sphere into concentric thin spherical shells, 
of thickness c?c, the potential of any one of these shells, of 
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radius c, at an internal point distant r from the centre will 

be obtained by writing c for 6, ^ for C, 47rcWc for Z7, in 

the first result of Art. 6. This gives 

p ^Tr<?do PfiY 47r P -p i j 

To obtain the potential of the whole shell, we must inte- 
grate this expression, with respect to c, between the limits 
a and a. This gives 

Again, the potential of the shell of radius c, at an external 
point, will be 






Integrating as before, we obtain for the potential of the 
whole shell, 

47r p p (a"^°-a'^') 

(2i + l)(2i + 3)A' * r*""' 

Suppose now that we wish to obtain the potential of the 
whole shell at a point forming a part of its mass, distant r 
from the centre. We shall obtain this by considering sepa- 
rately the two shells into which it may be divided, the 
external radius of the one, and the internal radius of the 
other, being each r. Writing r for a', in the first of the fore- 
going results, we obtain 

And writing r for a in the other result, we obtain 

(2i + l)(2i + 3) A* r**"- ' 

Adding these, we get for the potential of the whole 
sphere 

2i+ 1 A' t 2 *" "*■ (2i + 3) r*«J ' 
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It is hardly necessary to observe that the corresponding 
results for a solid sphere may be obtained from the foregoing, 
by putting a = 0. 

If the density, instead of being TiPiC*, be ^P^c^^ similar 

reasoning will give us, for the potential of the thin shell of 
radius c and thickness dc at an internal and external point 
respectively, 

And, integrating as before, we obtain for the potential of 
the whole shell, 

7 o » . IN/ ^-7-S\^-Pi(»"*"*'*^' — »'*""**') ^* at an internal point, 

47r p a"*"""*"^ — a *"'""^' 
(2i + l)(m + i + 3)/r^* ?"-' ^* ^ ^""^^"^^ P'"'''^- 

And, at a point forming a part of the mass, 



2i + l /r V m-i + 2 



„m-H+8 'm-H+5 -I \ 

^*'*" m + i + 3 P^j* 



11. Suppose, for example, that we wish to determine, in 
each of the three cases, the potential of a spherical shell 
whose external and internal radii are a, a', respectively, and 
whose density varies as the square of the distance from a 
diametral plane. 

Taking this plane as that of xy, the density may be ex- 
pressed by pz^, or pcV- Now fi^ = — ^ — . Hence the 
density of this sphere may be expressed as 

3A' » ■*"3 A*^ 

The several potentials due to the former term will be, 

2 

writing 2 for i and multiplying by « > 
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And for the latter term, writing for i, and 2 for m, and 
multiplying by g , 

47r p , , ,,. *!f A?LZiL" i'^ £. / Q^*~^* /-g'^^ 
12 A^^"* "^ ^' 15 K' r ' 3 A' V 4 "^ 5r > 

And, since P.r^ = — ^ — , we get for the potential at an 
internal point 

f. Iff (a' - O (3^' -r') + l («*-«")} ; 
at an external point 

at a point forming a part of the mass 

p f47r /g'-- r^ r^-a-^ 47r /g^^rV r^^g-^M 

A^IS I 2 + 7r' y^'^'^ "^^^ 3 I 4 + 6r /J ' 

12. We may now prove that by means of an infinite series 
of zonal harmonics we may express any function of fi what- 
ever, even a discontinuous function. Suppose, for instance, 
that we. wish to express a function which shall be equal to 
A from /Lt = l to /x = X, and to B from //»=X to /*== — 1. 
Consider what will be the potential of a spherical shell, 
radius c, of uniform thickness, whose density is equal to A 
for the part corresponding to values of fi between 1 and X, 
and to B for the part corresponding to values of fi between \ 
and — 1. 

Divide the shell, as before, into indefinitely narrow strips 
by parallel planes, the distance between any two successive 
planes being cdfi. 
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We have then, for the potential of such a sphere at any 
point of the axis^ distant z from the centre, 

for the first part of the sphere 



and for the latter part 

27rJ5c»Sc 



r 

•'-1 



dfi 



These are respectively equal to 
27r^c'Sc 



5£(p,+P,^ + P,^ + ... + P,^l+....)(Z/i, 



27rPc*Sc f^ 



£(p, + P,^ + P.J+... + P,?;+....)rf/., 



at an internal point ; and to 



z 



j\Po + P^l+-'+PiJi+-)dfM, 



at an external point. 

Now it follows from Chap. li. (Art. 23) that if i be any 
positive integer, 

whence, since I P,d/Lt = 0, it follows that 

•'-1 
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Also j P,dfJL=^l-\ j P,dfi=^l+\. 



Hence the above expressions severally become : 
For the potential at an internal point on the axis 

c 



[^(l-\) + 5(l+X)-^Y^{P.(X)-P.(X)} 



A-B,^,.. ^,...z- 



and for the potential at an external point on the axis 

Hence the potentials at a point situated anywhere are 
respectively 



27rc»& 



[{^(1-X)+5(1+\)}P» 



^-5 



{p,(x)-p,(\)}5Mir 



3 l-SVV -OVVJ g 



at an internal point; 
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and 

2ir<?hc \[A (1 - X) + J?(l + X)} -^^ 



^-^{P,(x)-P.(x)}^ 



at an external point. 

Now, if we inquire what will be the potential for the 
following distribution of density, 

i[A{l - X) + JS(l + X) - iA-B)[P^(\) - P,{\)}PM 

- (^-£){P.(x)-PXx)}P, (,*)-. .. 

- {A-B){P^,{X) -P,.,(X)}P.O*) - ...], 

we see by Art. 6 that it will be exactly the same, both at 
an internal and for an external point, as that above in- 
vestigated for the shell made up of two parts, whose densities 
are A and B respectively. * 

But it is known that there is one, and only one, dis- 
tribution of attracting matter over a given surface, which 
will produce a specified potential at every point, both ex- 
ternal and internal. Hence the above expression must 
represent exactly the same distribution of density. That is, 
writing the above series in a slightly different form, the 
expression 

^-^[x+ {P,(x) - Po(X)}P» 

+ {P.(x)-p,(x))P,o.*) + ... 

... + {P,„(X)-i^..(X)}P,(/t) + ...] 
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is equal to -4, for all values of /Lt from 1 to X, and to B for all 
values of [i from \ to — 1, 

13. By a similar process, any other discontinuous function, 
whose values are giyen for all values of [i from 1 to — 1, may 
be expressed. Suppose, for instance, we wish to express a 
function which is equal to A from //, = 1 'to ft = \, to jB from 
/L6 = \ to /A =s \, and to G from /lc = \ to /a = — 1, This will 
be obtained by adding the two series 

+ {P,j\)-p^,(\)}P,M + ...], 
+ {PhA)-p*.iW}p*W + •••]. 

For the former is equal to ^4 — jB from /lc = 1 to /^ = \, 
and to from /L6 = \to/Lt = — 1; and the latter is equal to 
B from /t* t= 1 to /L6 = Xj, and to G from /lc = \ to /lc =— 1. 

By supposing A and C each = 0, and jB = 1, we deduce a 
series which is equal to 1 for all values of fi from /a = \ to 
^ s= \^ and zero for all other values. This will be 



This may be verifiod by direct investigation of the 
potential of the portion of a homogeneous spherical shell, 
of density unity, comprised between two parallel planes, 
distant respectively c\ and c\^ from the centre of the 
spherical shell. 

14. In the case in which X^ and X, are indefinitely neaily 
eq[ual to each other, letXg = X, and\=X + d!X. We then 
have^ ultimately, 



i'.w-p.w-^^ 



F. a. 
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Hence P^,(\J - P^/\J - P^X,) ^ P^^(\J 

Hence the series 

^ {1 + 3P,(X)P» + 5P.(\)P» + ... 

+ (2i + l)P,(\)P,(^) + ...} 

is equal to 1 when fi = \ (or, more strictly, when fi has any 
value from \ to X + dX) and is equal to for all other values 

of /Lt. 

We hence infer that 

1 + 3P,{\)P,{fi) + . . . + (2; + l)P,{X)P,Oi) + ... 
is infinite when /a = \, and zero for all other values of fi. 

15. Representing the series 

i{l + 3P,{\)P,(fi) + ••. + (2*4- 1)P,(X)P^) + ...} 

by ij>{X) for the moment, we see that p(f>{\)d\ is equal to p 
when /t* = X, and to zero for all other values. Hence the 
expression 

{Pi<^(\)+P8<^(\) + ...]^ 

is equal to p^ when fi = \, to p, when fi == \,... Supposing 
now that \, \... are a series of values varying continuously 
from 1 to — 1, we see that this expression becomes 

I p(f>{\)dX, 

p being any function of \, continuous or discontinuous. 
Hence, writing ^(X) at length, we see that 

l[f_pdk + 3P» f pP,{\)d\ + ... 

+ (2{ + l)P,(ji)j.pPl\)dk + ...l 

is equal, for all ralues of fi from — 1 to + 1, to the same 
function of fi that p is of X. 






I* 
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16. The same conclusion may be arrived at as follows : 

< 

The potential of a spherical shell, "whose density is p, 
and volume Uj at any point on the axis of z, is 

pdK 
which is eqttal to g-j I prf\+-[ pPi(X)c?\+ ... 

for an internal pointy 

and to 21"" / P^ "*■ ~2 1 P-PiW^^ + ••* 

l| pP\{\)dK^.X 

for an external point. 

It hence follows that the potential, at a point situated 
anywhere, is 

for an internal point, '' 

and to ^IJ f ' pd\ + ^^f pPSK)dX + . . . 

for an external point. 

And these expressions are respectively equal to those 
for the potentials, at an internal and external point re- 
spectively, for matter distributed according to the following 
law of density : 
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+ .., 



+ (2i + l)P,(ji)j pPt(k)dX +...\. 

It will be observed, in applying this formula^ that if p be 
a discontinuous function of X, each of the expressions of the 

form I pPj[}C)dK will be the sum of the results of a series of 

integitU, each integration being taken through a series of 
values of \ for which p varies continuously. 



CHAPTER IV. 

SPHERICAL HARMONICS IN GENERAL. TESSERAL AND SEC- 
TORIAL HARMONICS. ZONAL HARMONICS WJTR THEIR 
AXIS IN ANY POSITION. POTENTIAL OP A SOUP NEARLY 
SPHERICAL IN FORM. 

1, We have hitherto discussed those solutions of the 
equation V*F= which are symmetrical about the axis of z, 
or in other words, those solutions of the equivalent equation in 
polar co-ordinates which are independent of ^. We propose, 
in the present Chapter, to consider the forms of spherical 
harmonics in general, understanding by a Solid Spherical 
Harmonic of the i^ degree a rational integral homogeneous 
function of x, y, z, of the i^ degree which satisfies the equa- 
tion V*F= 0, and by a Surface Spherical Harmonic of the 
1* degree the quotient obtained by dividing a Solid Sphe- 

rical Harmonic by («" +^ + «*)*• Such an expression, as we 
see by writing a? = r sin ^ cos ^, y = r sin ^ sin ^, « = r cos 0, 
will be of the i^ degree in sin 6 cos <^, sin sin 0, cos 6; and 
will satisfy the differential equation in Y^, 

or, "writing /a for cos ff, 

It will be convenient, before proceeding to investigate the 
algebraical forms of these expressions, to discuss some of 
their simpler physical properties. 

2, We will then proceed to shew how spherical har- 
monics may be employed to determine the potential, and 
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consequently the attraction, of a spherical shell of indefinitely 
small thickness. 

We will first establish an important theorem, connecting 
the potential of such a shell on an external point with that 
on a corresponding internal point. The theorem is as follows : 

If he the centre of such a shell, c its radim, P any irir 

temat pomt, P' an external point, so situated that P' lies on 

OP produced, and Hiat OP . OF = c", and if OP = r, OF = r', 

then the potential of the shell at P is to its potential at F 

* as c tor, or {which is the same thing) as r' to c. 

For, let A be the point where OP* meets the surface of 
the sphere, Q any other point of its surface. Then, by a 
known geometrical theorem, 

QP: QP iiAPiAFv.C'-r: f^^e. 

And 



c- 


— r 


cr- 


-f^ 


cr- 


or 







/ 


rr - 


-cr 


r'' 



Again, considering the element of the shell in the im- 
mediate neighbourhood of Q, its potential at P is to its 
potential at P' ^ QP is to QP, that is, as c to r, or (which 
is the same thing) as r' to c, which ratio, being independent 
of the position of Q, must be true for every element of the 
spherical shell, and therefore for the whole shell. Heace 
the proposition is proved. 

3. Now, suppose the law of density of the shell to be 

r 
such that its potential at any internal point is J?'(/Lt, ^) -^. 

c 

r* 
Then F(jm, ij>) -^ must be a solid harmonic of the degree t. 



Hence F(ji, ^) must be a surface harmonic of the degree e. 
Let us represent it by Y^. 

By the proposition just proved, the potential at any 
external point, distant r from the centre, must be 

Y -— 
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Hence, the compoiient of the attraction of the sphere on 
the internal point measured in the direction from the point 
inw.anlB,i.e, towards the centre of the sphere, is 

And the component in the same direction of the attraction 
on the external point, measured inwards, ia 

(*+i)y;^i. 

Now suppose the two points to lie on the same line 
passing through the centre of the sphere, and to be both 
indefinitely close to the surface of the sphere, so that r and r' 
are each indefinitely nearly equal to c. 

And the attraction on the external point exceeds the 
attraction on the internal point by 

(2i + l)-'. 
^ ' c 

Now, supposing the shell to be divided into two parts, 
by a plane passing through the internal point perpendicular 
to the line joining it with the centre, we see that the at- 
traction of the larger part of the shell on the two points will 
be ultimately the same, while the component attractions of 
the smaller portions, in the direction above considered, will 
be equal in magnitude and opposite in direction. Hence the 

Y 

difierence between these components, viz. (2i + 1) — , will be 

c 

equal to twice the component attraction of the smaller 

portion in the direction of the line joining the two points. 

But if Pi be the density of the shell, 8c its thickness, this 

component attraction is ivp^Bo. 

Y 
Hence (2i + 1) ^ = i^^p^ Sc, 

c 

2t + l „ 
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And, if cr^ be the correspondiDg surface density, 

21 + 1 



<r.= 



r. 



* 4c "'• 

It hence follows that if the potential of a spherical shell, 
of indejmitely small thickness^ he a surface harmonic, its 
potential cct any internal point will he proportional to the 
corresponding solid ham/ionic of positive degree, and its po- 
tential at any external point will he proportional to the 
corresponding solid harmonic of negative degree. 

That is, the propositioh* proved fot zonal harmonics in 
Chap. III. Art. 6, is now extended to spherical harmonics in 
general. 

4. The spherical hai^onic of the degree i will involve 
2i + 1 arbitrary constants. 

For the soHd spherical harmonic, r*F^, being a rational 
integral function of x, y, z of the i*** degree, will consist of 

-^ terms. Now the -expression V*F, being a 

rational integral function of x, y, z of the degree i — 2, will 

consist of -^ — ■—— terms ; and the condition that it must be 

• - (i—Vsi 

= for all values of x, y, z, will give rise to ^ — ^ — relations 

among the ^^ ^-^ coeflScients of these terms, leaving 

^^ ~ — -^ — 9~' ^^ ^** ^ ^> independent coefficients. 

5. We proceed to shew how the sphieridal harmonic of the 

degree i may be arranged in a series of terms, each of which 

may be deduced by differentiation' fifem the Zonal Harmonic 

symmetrical about the axis of z. Th^ solid zonal harmonic, 

which, in accordance with the notation already employed, is 

represented by r^P^ (ji), is a function of z and r of the degree i, 

d^V d^V d^V 
satisfying the equation V*F= 0, or --r-j + -r-,- 4- -tt = ©• 

Now, if we denote this lekpressioh by P^ (z), we see that 
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since it is a function of z and r, it is a function of the dis- 
tance {£) from a certain plane passing through the origin, and 
of the distance (r) from the origin. Further, if we write for z 
the distance from any other plane passing through the origin, 

dW cPV d^V 
leaving r unaltered, the equation ^-j + j^ + ~7^ =^ will 

continue to be satisfied*. 

Now z-V,OL (x + J—ly); a being any quantity whatever, 
represents the distance from a certain jplane passing through 
the origin, since in this expression, the sum of the squares 
of the coeflScients of z, x, y is equal to tmity. Hence 

P<{« + a(« + N/— ly)} is a solid zonal harmonic of the 

so 'U 

degree i, its axis being the imaginary line - = *J — = «. 
Therefore the equation 

da?^ dy"^ dz\ ^' 



is satisfied by F= P^ {« + a (a; + V— ly)}, that is, expanding 
by Taylor's Theorem, it is satisfied' by 

■*" 1.2... i dz' ' 

for all values of a. 

Hence, since the equation in V is linear, it follows that 
it is satisfied by each term separate^, 6t that, besides P^ {z) 
itself, each of the i expressions> 

satisfies the equation Y^Q. 

By similar reasoning we may shew that each of the i ex- 
pressions, 

satisfies the same equation. 
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Now each of the 2i solutions, thus obtained, is imaginary. 
But the sum of any two or more of them, or the result 
obtained by multiplying any two or more by any arbitrary 
quantities, and adding the results together, will also be a 
solution of the equation. Hence, adding each term of the 
first series to the corresponding term of the second, we A- 
tain a series of i real solutions of tha ^qoation. Another 
^uch iiarks nwy be obtemed by BuWacting each term of the 
second series fi-om the corresponding term of the first, and 

dividing by V— 1. We have thus obtained (including the 
original term Pi{z)) a. series of 2t + l independent solutions 
of the given equation, which will be the 2i + 1 independent 
solid harmonics of the degree i, 

6. We may deduce the surface harmonics from these by 
writing r sin cos for a?, r sin 6 sin ^ for y, r cos for z, 
and dividing by r*. Then, putting cos^ = //», and observing 

that P,(;5) = r*P,(/t), ^^ = r^^^^... we obtain the fol- 
lowing series of 2i + 1 solutions : 

cos<^sin5^:^, cos2^sin«^'^*l'^, .,, cosz>sin*5^^, 

sin^sin^^l^, 8in2<^sin»^^'^^, ... sin^<^sin'^^^^ 

Expressions of the form 

G cos or^ &m^d , \; , 

or Ssin a<j> &m''0 /^ , 

or their equivalents, 
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{0 and S denoting any quantities independent of 6 and <^) 
are called Tesseral Surface Harmonics of the degree i and 
order <t. The particular forms assumed by them when 
<r sf |R» called Sectorial Surface Harmonics of the degree i. 

It will be observed that, suoe — —-r-^ is a numerical constant, 

Sectorial Harmonics only involve 6 in Hie iaaoL 

8in*^,or(l-/x")*. 

The product obtained by multiplying a Tesseral or 
Sectorial Surface Harmonic of the degree i by r* (that is, 
the expression directly obtained in Art. 5) is called a Tesseral 
or Sectorial Solid Harmonic of the degree i. 

7, We shall denote the factor of a Tesseral or Sectorial 
Harmonic which is a function of 0, that is sin*'^ — r^-^ , or 

(1 —/**)' yl f ^J the symbol 7|f*^>, or, when it is necessary 

to particularize the quantity of which it is a function, by 
r,<'>0tA)or2;(<')(cosd). 

It will be ppnvenient, for the purpose of comparison with 
the forms of Tesseral Harmonics given in the Mecanique 
CHeste, and elsewhere, to obtain Tf''^ in a completely de- 
veloped form. 

Now, sinee P^iji) = «rT"«— o — • — ^i — % ^^ see that 
dy^' 2'. 1.2, 3.. .t c?/***" 

» 2i (2t -!)...(« -<r+l) /It*-' 
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- J (2i - 2) (2t -3)...(t -<T - 1) ^< 



-<r-i 



+ ^^-^ (2i - 4) (2» - 5)...(i - o— S)/*'-'-* 



•»•••» 



2»(2t-i)...(;-ff+i) |;^-'- ^iz|l^=Jp:Il/ti-.r-9 

(^♦- g-) (»• - ^ - 1) (^- - ^ - 2) (,'- ^ - 3) I 

■*" 2.4.(2i-l)(2t-3) '^ *■•;• 



And therefore 



''' - 2M.2.3,..i ^^ '^^ t'* 2(2»-l) '^ 

■*■ 2.4(2i-l){2t-3) '^ "•■•)• 

The form given by Laplace for a Tesseral Surface Har- 
monic of the de^ee i and order <r is (see Mecanique Celeste, 
Liv. 3, Cihapk 2, pp; 40—47) 

^ (1 - ,.r {/**-' - ^'^'lli^Zlf^^ f^-'-' + — } <=«« '^'^' 

-4 being a quantity independent of and <^ Tlie factor of 
this, involving ytt, is denoted by Thomson and Tait {Natural 
Philosopky, VoL 1, p. 149) by the symbol 0/*^>. Thomson 
and Tait also employ a symbol ^^^>, adopted by Maxwell in 
his Treatise on Electricity and Magnetism, Vol. 1, p. 164, 
which is equal to 



(^^-<7)(a^-o•-l),..(^-.c^ + l)^ ■ '^^ d/^^ 
'^(i + (r)(i+(r-l)...(t^-(r + l) * ' 



0*" - 1)^ {f^ 
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Heine represents the expression 

2 (2»- 1) *^ 

(^ - a) {i-a -l)(i-a- 2) (j - <r - 3) ._,., \ 
2.4.(2»-l)(2»-3) '^ ~-J» 

cr 

or (- 1)* 0/*^), by the symbol P^ (jjl), and calls these expres- 
sions by the name Zugeordnete Functionen Erster Art {Hand- 
luch der Kugelfwactionen, pp. 117, 118) which Todhunter 
translates by the term "Associated Functions of the First 
Kind," which we shall adopt, 

Heine also represents the series 

^ 2(2i-l) ^ 

(^-o•)(^-^-l)(t-^c^-2)(^-(^^^3) ,_ , 
■^ 2.4(2i-l)(2;-3) ^ 

by the symbol Ij^^fji), (p. 117). 

The several expressions, 3T['), 0|'), ^('), Pi, ^i, are con- 
nected together as follows : 

2M.2.3...t yj^j __ Q(^j 
2i(2i-l)...(i-(rH-l) * * 

- (.Ci'ir(<ti't)?.2.- M-'= (-!)'«- (1--VW- 

8. It has been already remarked that the roots of the 

equation P, = are all real. It follows also that those of the 

dP d^P 

equations -r-^— 0> -y-a == O.,, sxe real also. Hence we may 

arrive at the following conclusions, concerning the curves, 
traced on a sphere, which result from our putting any one 
of these series of spherical harmonics = 0. 

By putting a zonal harmonic=0, we obtain e small circles, 
whose planes are parallel to one another, perpendicular to 
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the axis of the zonal harmonic, and symmetrically sitiiated 
with respect to the diametral plane, perpendicular to this 
axis. If i be an odd number this diametral plane itself 
becomes one of the series. 

By putting the tesseral harmonic of the order cr=0, we 
obtain i — a small circles, situated as before, and <r great 
circles, determined by the equation costr^ = 0, or sin a^ = 0, 
as the case may be, iLeir plaiies all intersecting in the axis 
of the jsystem of liarmonics» the ac^e between the planes of 

any two tjonsecutive great circles being - . 

By putting the sectorial harmonic ^ 0, we obtain i 
great circles, whose planes all intersect in the axis of the 
system, the angle between any two consecutive planes being 

IT 






9. The tesseral harmonic may be regarded from another 
point of view. Suppose it is required to determine a solid 
harmonic of the degree t, and of the foim Yy, such that Y^ 
shall be the product of a function of fi, and of a function of ^, 
which ftmctioBS we will denote by the symbols if|, 4>|, respec- 
tively. The differential equation, to which this will Lead, is 

Now this will be satisfied, if we make J^ and 4>, satisfy 
the following two equations : 

The latter equation gives 

<I>^ = Ccos <r^ + (7 Bin o-^. 

And, taking o* aa an integer, positive or negativei the 
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former is satisfied by M,=^ T^^^\ i.e. (l^/^Wx^Y^'Cl -/*')'. 
as we proceed to prove. 

"We know that 

Differentiate a times, and we get 
whence, by Leibnitz's Theorem, 

or 

<r 

and, multiplying by (1 —/a')*, 

. +(.--<r)(t + <r+l)(l-.MT^'-0...(l). 

Now, putting (1-M*fp=2;<'). 
we get 
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rf/A* 



And «(i+l) 2'/') = i(»+l)(l-//f |5. 

+ (t-<r)(t + <r+l)(l-M')^p+<r'(l-/.')^-^' 

Hence the equation above given for M^ is satisfied by 
M^ = Ti^\ and the equation in T^ is satisfied by 

r,= or/') cos (r</> + (7T,<'>siti <7</>. 

10. In Chap. IL Art.lO we have established the fundamental 
property of Zonal Harmonics, that if i and m be two unequal 

positive integers, I P^PJI^l = 0, This is a particular case 

of the general theorem that if YJ, I^ be two surface har- 
monics of the degrees i and m respectively. 



flJ^T.YJi^^O, 
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For, let Vi, V^he the corresponding solid harmonics, so 
that Vt=r*Yt, F^ = r"'Y'„. Then, by the fundamental pro- 
perty of potential functions, we have at every point at which 
no attracting matter is situated, 

and therefore 

^'\da? ^ df ^ dz^) ^"^Kda? ^ df^ dz^)"^' 

or, in accordance with our notation, l^v"^«"" ^mV"'^ — 0. 

Now, integrate this expression throughout the whole 
space comprised within a sphere whose centre is the origin 
and radius a, a being so chosen that this sphere contains no 
attracting matter. We then have 

\\\{V,VV„-V^V'y)dxdydz = 0. 

But also, when the integration extends over all space 
comprised within any closed surfax^e, we have 

///(F.VT„- F.v'P^ ^dydz =//(f. ^ - vj^) dS^O^ 

dS denoting an element of the bounding surface, and -r- 
diflferentiation in the direction of the normal at any point. 

Now, in the present case, the bounding surface being a 
sphere of radius a, and V], V^ homogeneous functions of the 
degrees i w, respectively, 

d8^a'df.d<f>, ^^iar'T,. ^^^m-^F^, 

and, the integration being extended all over the surface of 
the sphere, the limits of ft are — 1 and 1, those of ^, and 2ir. 
Hence 

P.H. 6 
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•whence, ifm.-- i he not =0, 

The value of / 1 Y^d^id^ will be investigated here- 
after. 

11. We may hence prove that if a function of fi avd ^ 
can he developed in a series of mrfom harvfionics, stick de- 
velopment is possible in only one way. 

For suppose, if possible, that there are two such develop- 
ments, so that 

and also 

Then subtracting, we have 

o=r^-F;+r,-r/+... + r,-r;' + ... identically. 

Now, each of the expressions Y^ — T/, J^ — F/... rj — F/ 
being the difference of two surface harmonics of the degree 
0, 1, •••i... is itself a surface harmonic of the degree 
0, 1, ...t.... Denote these expressions for shortness by 
^, Zj...-^... so that 

=^0 + ^ + — + 4+ — i^ei^ically. 

Then, multiplying by Z^ and integrating all over the 
surface of the sphere, we have 



0=r r'z^dpd^. 



That is, the sum of an infinitQ number of essentially 
positive quantities is =0. This can only take place when 
each of the quantities is separately == 0. Hence ^, is identi- 
cally = 0, or Y^ ~ Fp and therefore the two developments 
are identical 

We have not assumed here that such a development is 
always possible. That it is so, will be shewn hereafter. 
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12. By referring to the expression for a surface har- 
monic given in Art, 4, we see that each of the Tesseral and 

Sectorial Harmonics involves (1 — ft*)*, or some power of 

(1 — ft*)*, as a factor, and therefore is equal to when /Lt = ± 1. 
From this it follows that when /it = ± 1, the value of the 
Surface Harmonic is independent of ^, or that if Y {/m, (p) repre- 
sent a general surface harmonic, Y{± 1, ^) is independent of 
<l>, and may therefore be written as F (+ 1). Or F(l) is the 
value of xQi, if>) at the pole of the zonal harmonic PiQi), 
ir(— 1) at the other extremity of the axis of P, (ji)^ 

We may now prove that 



r.#=2,rr,(i)P.(/^). 



I 

For, recurring to the fundamental equation^ 

Now, if we integrate this equation with respect to if>, 
between the limits and 27r, we see that, since 

and the value of rj only involves (j) under the form of cosines 
or sines of d> and its multiples, and therefore the values of 

-^ are the same at both limits, it follows that 
Hence 

- . |i{(^-'*')(/o"w)}+i(i+i)(f;w)=o. 

Hence | Y^d(f> is a function of fi which satisfies the 

fundamental equation for a zonal harmonic^ and we therefore 
have 

6—2 
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r 





C being a constant, as yet unknown. 

To determine O, put ft=l, then by the remark just made, 
Yf becomes 3^4(1), and is independent of (f). Hence, when 

/[* = !,[ r,(Z^=:27rr,(l). Also P,(/a) = 1. We have there- 
fore ^ 27rr,(l) = a, 

Jo 
It follows from this that 

13. "We may now enquire what will be the value of 

ri r2n 
j J Y,Z,dfid4>, 

Y^y Z^ being two general surface harmonics of the degree i. 
Suppose each to be arranged in a series consisting of the 
zonal harmonic P^ whose axis is the axis of z, and the system 
of tesseral and sectorial harmonics deduced from it. Let us 
represent them as follows : 

+ C^Tp cos ^ + C,!r/«) cos 2^ + ... + C^Ty^ cos o-^ + ... 

+ 0,Tp cos ul> 
+ S^TPsm(l> + SJP 8m2<l> + ... + S^T^^^^ sma<f> + ... 

Z,=^ aP, 

+ c^Tp cos ^ + c,r/2) cos 2(f> + ... + c^r/<^) cos <r^ + ..• 

C^T^^''^ COS %<!> 

+ 8^Tpsm4> + 8jrj^ sin 2^ + ... +8^T^^^^ sm<7^+ ... 

Hence the product Y^Z^ will consist of a series of terms, 
in which ^ will enter under the form cos tr^ cos a^, or 
COB (T<j> sin a' (f). This expression when integrated between 
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the limits and 27r vanishes in all cases, except when 
<r' = <r and the expression consequently becomes equal to 
cos' (r<}), or sin\(r^. In these cases we know that, a being any 
positive integer, 

I cos' <7^ d^ == I sin' <r(f>d^ = tt. 

JO •'0 

Hence the question is reduced to the determination of the 
value of 

Now ?;w=(i-^«)^^ • 

~2M.2.3...i^ '^^ dfi**' ' 

But, by the theorem of Rodrigues, proved in Chap. II. 
Art 8, we know that 

Hence T^^'^ may also be expressed under the form 

whence it follows that 

^ ' '~^ ^ 12' J. 2.3.. .i; i-a- d/*<+' dfi*-' ' 



Now, putting {/I' — lY'siM for the moment, and inte- 
grating by parts, 

f d^+'Md*-'M , _ d'+'-^ M d ^-' M 
) dfi*+' d/ju*-' '*" dft**'-^ ~djly^ 



-/ 
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The factor -j—. — vanishes at both limits, hence 
dfi*'" 

fi di+'Md*-'M j __ fi df+'-^M di-'+-^ M , 

- (-i\* P di*'-' M d^-'+» M J 

by a repetition of the same process. 

And by repeating this process 0- times, we see that 
fi. d**'Md*-'M, , -,, p /(?il/V , 

= (- 1)' (2'. 1 . 2 . 3 ...i7 r P,*rf/* 
= (-.l)'(2M.2.3...ir2-|Y. 

and therefore 

r f*'^/'>oos<r^)*(iAMZ^=f /^"(T.Wsino-^/rf/id^ 

It will be observed that this result does not hold when 
o- = 0^ in which case we have 

Hence I / Y^Z^dyid^ 



/I I i + <r 9 



2i + l 



* In this ease f toa*vfd^='f 8in*«'^({^=2r. 
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" 27r fU' + l |*+2 

14 We have hitherto considered the Zonal Harmonic 
under its simplest form, that of a " Legendre's Coefficient " in 
which the axis of z, i.e. the line from which is measured, is 
the axis of the system. We shall now proceed to consider it 
under the more general form of a "Laplace's Coefficient," 
in which the axis of the system of zonal harmonics is in any 
position whatever, and shall shew how this general form may 
be expressed in terms of Ft{fi) and of the system of Tesseral 
and Sectorial Harmonics deduced from it. 

Suppose that ^, ^' are the angular co-ordinates of the 
axis of the Zonal Harmonic, i.e. that the angle between this 
axis and the axis of z is ff^, and that the plane containing 
these two axes is inclined to a fixed plane through the axis 
of z which we may consider as that of zx, at the angle ^\ 
In accordance with the notation already employed, we shall 
represent cos ff by fi. 

The rectangular equations of the axis of this system 
will be 

X _ y __ z 

sin 0" cos (j)' "" sin ^ sin <^' "" cos ^ ' 

Hence the Solid Zonal Harmonic of which thia is the axis 
is deduced from the ordinary form of the solid zonal bar- 
monic expressed as a function of z and r by writing, in plaoe 
o£z, xbwlO' cos<l>' + y sinff ain^' + zcosff. 

To deduce the Surface Zonal Harmonic, transform the solid 
zonal harmonic to polar co-ordinates, by writing r sin ^ cos ^ 
for 0?, r sin d sin for y, r cos for z, and divide by r*. 

The transformation from the special to the general 
form of surface zonal harmonic may be at once eflFected, 
by substituting for fi,OT cos^, cos^cosa+sin5sin^cos(^— ^'). 

Now, in order to develope 

Pf {cos ^ cos ^ 4- sin ^ sin ^ cos (^ — <f>)] 
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in the manner already pointed out, assume 

P^ {cos ^ COS ^ + sin sin 0" cos {<}) — <l>)} 

= ^P,0i) + (O(i)cos^ + ^Wsin</>) Tp 

+ (C(2) cos 2<j> + iS(2) sin 2^) 2^(2) + ... 

+ ( (?<*') cos o-^ + iS<<^> sin o-^) r/*') +.. . 
+ ( (7<»^ cos i<l> + 8^^ sin i<f>) T(f>, 

the letters j4, ... (7^*^^ iS<')... denoting functions of fi and 
i^', to be determined. 

To determine C^^^V niultiply both sides of this equation 
by cos a^Ty^ and integrate all over the surface of the sphere, 
ie. between the limits — 1 and 1 of fi, and and Zir of <j>. 
We then get 

I jPt {cos cos ^' + sin ^ sin ff cos (^ - if>)} cos o-^ Z;<*^> rffw?^ 

= (7 ^r J*' (cos (r^Z;<<^>)" dfid<l> 

It remains to find the value of the left-hand member of 
this equation. 

Now cos <7<^2y<'> is a surface harmonic of the degree i, and 
therefore a function of the kind denoted by Y^ in Art. 12. 

And we have shewn, in that Article, that 

4 

that is, that if any surface harmonic of the degree i be multi- 
plied hy the zonal harmonic of the same degree^ and the prodvxit 
integrated all over the surface of the sphere, the integral is 

equ>al to ^. — j into the value which the surface harmonic 

assuToes at the pole of the zonal harmonic. 
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Hence 

'1 raw 



I I P,{cos^cos^ + sm^sm^cos(^ — ^^}r, {fi, (j)) dfid<j> 

and therefore 

j f Pi {cos ^ cos ^ + sin sin 0' cos {(f> - ^')} cos a^T^'^dfidtf) 

Hence 

2z + l T i yr'j U-o-21 + 1 ' 



I— o- 



or C^*') = 2 h==cos <rA' Tjf^ (ji). 



Similarly ;SfW = 2 ^4=^ sin o-f 3;<<^> 0*')- 



And to determine ji, we have 
r pP, {cos ^ cos ^ + sin ^ sm ^ cos (^ - ^')} ^i W ^M# 

47r p . /\ _ J 47r 

or ^ = P,0*'). 
Hence, P, {cos ^ cos ^ + sin sin ^ cos {<j> - ^')} 

= p. OjO p, 0*) + 2 l^cos (^ - <!>•) T/i) 0*0 ?.« (m) 

+ 2 [=1 cos 2 (.^ - ^') T(f> 0*') 2;® 0*) + . . . 
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+ 2l|=cos<r (^-f ) Tyi iij!) T^'^iji) + ... 



+ 2 j^ COS t (* - f ) 2;» 0*') JiW 0*). 

15. We have already seen (Chap. li. Art. 20) how any 
rational integral function of fi can be expressed by a finite 
series of zonal harmonics. We shall now shew how any 
rational integral function of cos^, sin ^ cos ^, sin sin ^, 
can be expressed by a finite series of zonal, tesseral, and 
sectorial harmonics. 

For any power of cos j> or sin ^, or any product of sutsh 
powers, may be expressed as the sum of a series of terms of 
the form cos <r<f>, or sin o-^, the greatest value of a- being the 
sum of the indices of cos ^ and sin if>, and the other values 
diminishing by 2 in each successive term. Hence any 
rational integral function of cos 6, sin d cos ^, sin d sin ^, will 
consist of a series of terms of the form 

cos*" 6 sin* 6 cos aj> or cos'" 6 sin* 6 sin <r^, 

where n is not less than <r. 

If w be greater than c, w — <r must be an even integer. Let 
w — (7 = 25, then writing sin"^ under the form (1 — cos"^)' sin""^, 
we reduce cos"* sin* cos aj> to the sum of a series of terms 
of the form cos^ sin*' cos aj>, or, writing cos = fi,oi the 

form fi^ (1 — fi*) * cos a^. 

Similarly cos** ^ sin* ^ sin <r^ is reduced to a series of 

terms of the form fi^ (1 — ^') ' sin cr^. 

1 d^ 

and fiP"^^ can be developed in a series of terms of the form 
of multiples of Pp+«r, Pp+«r-2 ... . (Chap. IL Art. 17.) 

Hence fi^ can be expressed in a series of the form 
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Af^y A^ representing known numerical constants^ and therefore 
/aP (1 — fi^^ assumes the form 

consequently multiplying these series by cos a^ or sin <r^, we 
obtain the developments of 

— J? 

/ttP (1 — /tA*) ' coso*^ and /LtP(l — /tA')*sin<r^ 

in series of tesseral harmonics. 

16. We will give two illustrations of this transformation. 

First, suppose it is required to express cos' sin' sin ^ cos ^ 
in a series of Spherical Harmonics. 

Here we have sin ^ cos ^ = - sin 2^. 

Hence cos' sin' sin ^ cos ^ = ^ cos' sin' sin 2^. 

Comparing this with cos"* sin* sin <r^, we see that n is 
not greater than <r. 

Hence cos' sin' sin ^ cos ^ = ^ /Lt' (1 — /Lt') sin 2<f). 



^^ '^'"O^'** 



af> 



and . M* = J-P« + |-P, + |i'., 

••'^ 12 \35 d/t* ^ 7 d/tV 

105 rf/*' "^ 21 d/t" ' 
.'. cos* ^ sin* d sin ^ cos ^ 
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Next, let it be required to transform cos'l? sin' 5 sin ^ cos'<^ 
into a series of Spherical Harmonics, 

Here sin ^ cos* ^ = ^ sin 2^ cos ^ = t (^^^ ^^ + ^^^ ^)* 
Now cos' 6 sin* 6 sin 3^ = /tt* (1 - fj?)^ sin 3<^ 

Also cos' ^ sin* ^ sin i>=iJ? (1 - ij^) (1 — /t')isin^ 

= (/t»-/)(l-/t')isin^ 



= {i|'**-B|'*')(^-'^')'«^^- 



Also (Chap. n. Art. 17) 

'^ 35 * 7 • 5 •' 



Hence cos* 6 sin* ^ sin 3^ 

1 /16<?P, , 24<?PA,, ... . _. 
= 120(231^ +77^j(^-''^'^"^^^ 

Andcos*^sm'^sin^ = -^gg3-^«+.5^^^ + g3-^« 

2 dP. IrfP.N,, ^1 . . 

:s-(—TO)-~TO-) iydAsinit 

V693 • 385 * 63 » J'^^"-'P> 

.-. coa*^sin*^sin^cos*^=|g^Z;(») + j^2';«}sin3^ 

_ {A. yd). J-j- a)_JL TM sin*. 
1693 • 770 * 63 * J ^ 
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17. The process above investigated is probably the most 
convenient one when the object is to transform any finite 
algebraical function of cos 0, sin ^ cos <^, and sin sin ^, into 
a series of spherical harmonics. For general forms of a 
ifunction of fi and ^, however, this method is inapplicable, 
and we proceed to investigate a process which will apply 
imiversally, even if the function to be transformed be discon- 
tinuous. 

We must first discuss the following problem. 

To determine the potential of a spherical shell whose 
surface density is F(ja,^), JP^ denoting any function whatever 
of finite magnitude, at an external or internal point. 

Let c be the radius of the sphere, r' the distance of the 
point from its centre, ^, ^' its angular co-ordinates, V the 
potential. Then fi being equal to cos ^ 



J -1 ; [r'"- 2cr' {cos cos 0' + sin ^ sin 0' cos (<^ - 0')} + c"]* 

The denominator, when expanded in a series of general 
^onal harmonics, or Laplace's coefficients, becomes 

for an internal and an external point respectively, P, (/Lt, ^)- 
beiug written for 

P^ {cos ^ cos ^ + sin ^ sin ff cos (^ — ff>) }. 

Hence, F^ denoting the potential at an internal, F, at an 
external, point, 
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It will be observed that the expression P^ (ji, <f)) involves 
fi and fi symmetrically, and also ^ and if>\ Hence it satisfies 
the equation 

And, since fi and <h are independent of fi and ^', this 
differential equation will continue to be satisfied after P, has 
been multiplied by any function of fi and ^, and integrated 
with respect to fi and <f>. That is, every expression of the 
form 

fJl'PiO^^ 4!) F(ji, 4>) dii.d4> 

is a Spherical Surface Harmonic, or "Laplace's Function" 
with respect to ji and ^' of the degree i. And the several 
terms of the developments of Fj are solid harmonics of the 
degree 0, 1, 2,.,i... while those of F, are the corresponding 
functions of the degrees — 1, —2, — 3... — (i+ 1), ... And 
these are the expressions for the potential at a point (r', /tt', ff>) 
of the distribution of density F(fi\ <}>') at a point (c, /a, ((>). 

Now, the expressions for the potentials, both external 
and internal, given in the last Article, are precisely the same 
as those for the distribution of matter whose surface density is 

+ {2i + l)fJ^^P,0i,i>)F(fL,i>)dfi,di>+..j, 
or, as it may now be better expressed^ 
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+ 3 J [ P^[cos0cos0'+ sine smd" COS (ij>''4})FQi,4>)dfid4> 
+ ... 
(2i+l) [ j P,{cos^cos^+sin^sm^co9(^^')}F(/i,^)d/Arf</)+... . 

And, since there is only one distribution of density which 
will produce a given potential at every point both external 
and internal, it follows that this series must be identical 
with F(jj!, ^y We have thus, therefore, investigated the 
development of F(u\ d>) in a series of spherical surface 
harmonica* 

The only limitation on the generality of the function 
F(ji\ <f)) is that it should not become infinite for any pair of 
valuer comprised between the limits —1 and 1 of /a, and 
and 27r of ^. 

18. Ex. To express cos 2(f> in a series of spherical har- 
monics. 

For this purpose, it is necessary to determine the value of 
(2/+ 1)1 i P^ {cos^ cos^+sin^ sin^ cos (^-^')} ^^^ 2<f>d/id4>. 

Now P, {cos ^ cos ^ + sin sin ff cos (^ — (f>)} 

= P,(cos^)P,(cos^') 

^ 2 . ^dP,(eos^ . ^dP,(cos^) ., .,, 

2 

+ (i-l) »(» + !)(* + 2) ^""'^ 

<?P, (cos ^ . ,^ <?P, (cos ff) „ , . ■„ ^ 
Now I cos<r(^ — ^')co92^d^ = 0, 

JO 

for all values of cr except 2. 

* In connection with the subject of this Article, see a paper by lir G. H. 
Parwin in the Messenger of Mathematics for March, 1877. 
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And I cos2(^ — ^')cos2^(^aswco8 2^'. 
Jo 

Also 
And 

^1 (ifr*— 1)* 

s:2M • 2 . 3...i/iP| ^-i=i . 

Now when /tt = 1, 

And when /* = — !, 

Hence 
/\in«^gc?;. = ,,^j3 . 2M.2.3...qi-(-in 

= 4 or 0, as i is even or odd ; 

.-. r J*" sin'g ^-^JT^^ cos 2 (^ - f ) cos 2(^ dfid^ 

s 47r cos 2^' or 0, as » is even or odd ; 
.•. cos 2^' 

1 fc 2 , . ,„d»P,(co8^) -., 
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2 , . , ^ d'P, (cos 5') ^., 

4 13 .. ^ >, Q 4smV ^j-7i — ^7rcos2^ 

5.6.7.8 d/Lt' ^ 

^ } 

= 2 cos 2f (^, + 3^g + ^1|^3 + ....). 

Hence the potential of a spherical shell, of radius c and 
surface density cos 2^', will be 

and 

at an internal and external point respectively. 

19. We will now explain the application of Spherical 
Harmonics to the determination of the potential of a homo- 
geneous solid, nearly spherical in form. The following 
investigation is taken from the M4canique Celeste, Liv. iii. 
Chap. II. 

Let r be the radius vector of such a solid, and let 

r = a + a (a^ Fj + a, F, H- ... + a< Y] + ...), 

a being a small quantity, whose square and higher powers 
may be neglected, a^, aj,,...a,... lines of arbitrary length, and 
r^, F,,... F,... surface harmonics of the order 1, 2,...t... re- 
spectively. 

The volume of the solid will be « 7ra\ 

For it is equal to 

I r^drdfMd<l> 

1 /■' f2» 



4 , 



', since 1 1 Yidfid^ = 0, 



for all values of i. 
F. H. 
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Again, if the centre of gravity of the solid be taken as 
origin, a^ = 0. 

For if z be the distance of the centre of gravity from the 
plane of xi/, 



I 1 I 7^fidrd/jLd<f> 



-\'4ia^2 {a^l\^\•a^l\'^...'^c^Y^'\^ ...)dfid^ 





= 4«a'u a^\ I /A Y^dfid(l>, 



Similarly 
4 



3 """ 



X = 4a' a . a^ I I (1 - /x')* cos ^ F^ rf/xc?^, 

4 r^ r^' 

^ Tra'y = 4a' a . a^ I I (1 — fx*)i sin ^F^ dfid<f}. 

Now Fj is an expression of the form 

A/i + Bil" fi*)^ cos ^ + 0(1 - /jT)^ sin ^, 

and therefore all the expressions x, y, z cannot be equal to 0, 
unless a^ = 0. 

We may therefore, taking the centre of gravity as origin, 
write 

r=^a + a{aj^Y^'{'... 4- a^Yi+ ...), 

as the equation of the bounding surface of the solid. 

Now this solid may be considered as made up of a homo- 
geneous sphere, radius a, and of a shell, whose thickness is 

a(a,Fj4- ... + a,l^ + ...). 

The potential of this shell, at least at points whose least 
distance from it is considerable compared with its thickness, 
will be the same as that of a shell whose thickness is aa, and 
density 



^•(^'■»+- + ?^'+--)' 
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Pg being the density of the solid. Therefore the potential, 
for any external point, distant R from the centre, will be 

*m 3^ + 4m»« [-Y E' + "-+ ^i^in ^ + •• -J • 

The potential at any internal point, distant jB from the 
centre, will be made up of the two portions 

I irp.R' + 2irp, (a' - iJ*) or 2vp (a* - y) 

for the homogeneous sphere, 

for the shell, and will therefore be equal to 

20. If the solid, instead of being homogeneous, be made 
up of strata of different densities, the strata beiDg concentric, 
and similar to the bounding surface of the solid, we may 

deduce an expression for its potential as follows. Let - r be 

the radius vector of any stratum, p its density, r having the 
same value as in the last Article, and p being a function 
of c only. Then, he being the mean thickness of the stratum, 
that is the difference between the values of c for its inner 

« 

and outer surfaces, the potential of the stratum at an ex- 
ternal point will be 






To obtain the potential of the whole solid at an external 
point we must integrate this expression with respect to c, 
between the limits and a, remembering that p is a func- 
tion of c. 

7—^ 
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Again, the potential of the stratum, above considered, 
at an internal point will be 

47rpcSc + 4^pa — (^-Y«-^ + -^ -^ + ... 

To obtain the potential of the whole solid at an internal 
point we must integrate the expression (1) with respect to c 
between the limits and B, and the expression (2) with 
respect to c between the limits It and a, remembering in 
both cases that p is a function of c, and add the results 
together. 



CHAPTER V. 

SPHERICAL HARMONICS OF THE SECOND KIND. 

1. We have already seen (Chap. il. Art. 2) that the 
differential equation of which P^ is one solution, being of 
the second order, admits of another solution, viz, 

fip f ^^ 

Now if fi between the limits of integration be equal 
to ± 1, or to any roots of the equation P, = (all of which 
roots lie between 1 and — 1), the expression under the 
integral sign becomes infinite between the limits of inte- 
gration. We can therefore only assign an intelligible 
meaning to this integral, by supposing fi to be always be- 
tween 1 and 00 , or between — 1 and — oo . We will adopt 
the former supposition, and if we then put C7= — 1, the 

expression p,,^ _ ^. (i.e. p2 / a^iJ ^^^ ^® always posi- 
tive. We may therefore define the expression 



I. 



as the zonal harmonic of the second kind, which we shall 
denote by Q^y or Q^ (ji), when it is necessary to specify the 
variables of which it is a function. 

It will be observed that, if fi be greater than 1, P^ is 
always positive. Hence, on the same supposition, Q^ is 
always positive. 

Weseethat $,= | -r^ = 9^^g— T> 
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^ r da 

Aad, in a similar manner, the values of Q,, Q^,.>. may 
be calculated. 

2. But there is another manner of arriving at these 
functions, which will enable us to express them, when the 
variable is greater than unity, in a converging series, with- 
out the necessity of integration. 

This we shall do in the following manner. 

I^t £/"« , V being not less, and /a not greater, than 

V — fl 

unity. 

Then ^=^_1„ £^„_^, 

, dp ('' — /*)*' dfjk (p^fi)*' 

±ifi-^^-J^^fJ!L 2 \ ^ 1-fip 

dpX" ^ dvj~(v-fi')V-l v-fi)~ (v-fi)*' 

•••^(a-'4^-i{a-'")f}- 

1 . • 

Now, let be expanded in a series of zonal bannonics 

P», P»w.P,0i),8otbat 

Then|^|(l-M')^^=...-t(.-+l)^.(.')P» + ... 
by the definition of P (/a). 
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And these two expressions are equaJL Hence, equating 
the coefficients of P< (/i.), 

Hence <^,(^) satisfies the same diflferential equation as P, 
and Q,. But since f7= when i/ = oo , it follows that <j>i{v)=0 
when 1/ = 00 . Hence <^, {v) is some multiple of Q^ {v) =-4Q, (i/) 
suppose. It remains to determine A. 

Now, ^,(1^) may be developed in a series proceeding by 
ascending powers of - , as follows. 

We have = - + ^ + ...+-^ + -... 

V — fM V IT V 

and also = ^^W -Po(m) + ^iM ^M +- + ^< W ^*(m) + --' 

Now, by Chap. ii. Art. 17, we see that, if m be any 
integer greater than i, the coefficient of P^ in /a"* is 

(2^ + l) , <"^-:+^) ^"?~:^'^^^••^"lw\^^ if *be odd. 
^ ^ (m + t + l)(m + t-l)...(»» + 4)(m + 2) 

and (2^+1) , ^^•r,Tf''5^"ir'T*ln;r..i. ^^^^^^^'^' 

^ ^ (m+t+1) (m + « — l)...(wi + 3)(w+l) 

m — 1 being always even. 

Hence, writing for m successively i, i + 2, / + 4, .•• we get 

*.(") = I2t + 1) |(2i+l)(2i-l)...(i + 2) v^' 

4.6...(t + l) 1 



(2t + 3H2i + l)...(i + 4) v' 



.<+t 



+ (2 



2?T 5)(2/+3)...(t + 6) ,r-V+-j ^^^^ °'^<^' 
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and =(2»-+l){pj:p 



2. 4...t 



l)(2t-l),.. (» + !)„« 
4. 6.. .(1 + 2) 1 



<+8 



6.8...(; + 4) 1 



' (2i+5)(2t + 3)...(t + 5) J/' 
Now, recurring to the equation 



i+6 



if i be even. 






we see that, if Q<(i;) be developed in a series of ascending 

1 1 

powers of -, the first term will be ^ .,, . — tt-ht, where C 

is the coefficient of /jl* in the development of P^ (/i) ; 
^, . ^ (t+2)(t + 4) ...(2i-l) .... ,, 



and 



^ (t + lHi + 3) (t + 5)^,,(2i^-l) 

2.4.6...i^ "^^^^ 



even. 



Hence the first term in the development of Q^ {li} is 

2.4.6...(;-l) . 



and = TT 



(i + 2)(t + 4)...(2i-lj (2i + l) 
2.4.6...i 



(i+l)(i + 3)...(2i-l)(2i + l) 



if i be odd, 



if i be even. 



which is the same as the first term of the development of 
P, (i/), divided ^7 2i+i * 

Hence A = 2i+ 1, and we have 
1 



V — fA 



= Co W Po W +3(2, M P» + 5Q,W P» + ... 



3. The expression for Q^ may be thrown into a more 
convenient form, by introducing into the numerator and de- 
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nominator of the coefficient of each term, the factor neces- 
sary to make the numerator the product of i consecutive 
integers. We shall thus make the denominator the product 
of i consecutive odd integers, and may write 

^^^""^ " I .S . b..,{2i+l) v*^''^ S . 5 .7...\2i + 3) p*^ 






5.6.7...(t + 4) 1 
5.7.9...(2i+5)i/* 

(24 + 1) (2A?4-2)...(i + 2A?) 



,»!+ .... 



" {2fc+l) (2A; + 3)...(2i + 2A; + l) i/*^ 
whether i be odd or even. 

4. We shall not enter into a full discussion of the pro- 
perties of Zonal Harmonics of the Second Kind. They will be 
found very completely treated by Heine, in his Handbuch der 
Kugel/unctionen. We will however, as an example, investi- 

dQ,, 

dp 



gate the expression for -^ in terms of Q^^, ^^^.g... 



Recurring to the equation 

^= <?.(") P» + 3<2,(i') P» + . .. 

+ (2t + l) 5(2.(1-) P» + ... 
we see that 



dfi 



Now we have seen (Chap. ii. Art. 22) that 
^^ = (2» - 1) P... (,*) + (2t - 5) P... 0*) + ... 

Hence ^%^ = (2» + 1) P,(/*) + (2t- 3) P^iji) +... 
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^^^=(2i + 5)P^(^)+(2i + l)P» + .., 

^^ = (2;+ 9) P^M + (2i + 5) P^Ox) 

+ (2i+l)P,0) + ... 

And therefore the coeflScient of P^Ott) in the expansion 
p d 1 , 

01 -3 IS 

afjb v — fi 

(2i+l) {(2i+3) Q^,W + (2t+7) <?.«(!') + (2i+ll) Q^{v) + ...]. 
Again, 

+ (2i + l)^P/^). 

And T h J = 0. 

av V — ii> aiMv—/i 

Hence, comparing coefficients of P, (/t), 

^^^ = -(2i + 3)Q,,.M-(2t + 7) (?^M 

-(2t + ll)(2.„(r)-... 
Hence it follows that 

and therefore that 



/, 



"<?.„ W rf*- = 2?T3 l^^'W-Qu.W}. 



5. By similar reasoning to that by which the existence of 
Tesseral Harmonics was established, we may prove that there 
is a system of functions, which may be called Tesseral Har- 
monics of the Second Kind, derived from r/'> in the same 
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manner as Q, is derived from P^. The general type of such 
expressions will be 

Ti-)iv) r ^^ 

and this when multiplied by cos<T<^or sino-^, will give an 
expression satisfying the differential equation 

|(l-/*«)|^yZ7+{t(i + l)(l-^«)-a*}U'=0, 

and which may be called the Tesseral Harmonic of the 
second kind, of the degree i and order cr. 



CHAPTER VI. 



ELLIPSOIDAL AND SPHEBOIDAL HARMONICS. 



1. The characteristic property of Spherical Harmonics 
is thus stated by Thomson and Tait (p. 400, Art 537). 

"A spherical harmonic distribution of density on a spheri- 
cal surface produces a similar and similarly placed spherical 
harmonic distribution of potential over every concentric 
spherical surface through space, external and internal." 

The object of the present chapter is to establish the exr 
istence of certain functions which possess an analogous pro- 
perty for an ellipsoid. They have been treated of by Lam^, 
in his Lefons sur les fonctions inverses des transcendantes et 
les fonctions isothermes, and were virtually introduced by 
Green, in his memoir On the Determination of the Exterior 
and Interior Attractions of Ellipsoids of Variable DensitieSy 
(Transactions of the Cambridge Philosophical Society, 1835). 
We shall consider them both as functions of the elliptic co- 
ordinates (as Lam^ has done) and also as functions of the 
ordinary rectangular co-ordinates; and after investigating 
some of their more important general properties, shall pro- 
ceed to a more detailed discussion of the forms which they 
assume, when the ellipsoid is a surface of revolution. 

2. For this purpose, it will be necessary to transform 
the equation 

into its equivalent, when the elliptic co-ordinates e, v, v* are 
taken as independent variables. If a, h, c be the semiaxes 
of the ellipsoid, the two sets of independent variables are 
connected by the relations 
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Thus a' +6, J' + e, c' + e are the squares on the semiaxes 
of the confocal ellipsoid passing through the point a?, y, z, 

c? + V, 6' + V, (? + V, the squares on the semiaxes of the 
confocal hyperboloid of one sheet. 

a' + v\ V 4- v', c* + Vy the squares on the semiaxes of the 
confocal hyperboloid of two sheets. 

Thus, € is positive if the point a?, y, z be external to the 
given ellipsoid, negative if it be internal. 

And, if a' be the greatest, c' the least, of the quantities 
a\ h\ c«, 

6 will lie between — c* and oo , 

^ w a ^ n ^ f 

V „ „ -a* „ ^b\ 

3. Now -7-j + -i-i + Tn = is the condition that 
dor axf dzr 

taken throughout a certain region of space, should be a mini- 
mum. In the memoir by Green, above referred to, this 
expression is transformed into its equivalent in terms of a 
new system of independent variables, and the methods of the 
Calculus of Variations are then applied to make the resulting 
expression a minimum. We shall adopt a direct mode of 
transformation, as follows : 

Suppose a, )8, 7 to be three functions of x, y, z, such that 

V»a = 0, v'/3 = 0, vV = (1), 

such also that the three families of surfaces represented by 
the equations a = constant, fi = constant, 7 = constant, inter- 
sect each other everywhere at right angles, i.e. such that 
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dfidy dfidy d/3dy_ dy d^ dydx cJydi^^ 
dx dx dy dy dz dz ' dx dx dy dy dzdz~ ' 

d2d0^d2d0_^d2d0^^ 

dx dx dy dy dz dz ^*'^' 

Then 

dV^dVda dVd^ dV&y 
dx di dx d^ dx djy dx^ 

ds^'~U\dx) ^ d/a^ydxj "^ dy* \dx) 

^d^d^^ ^d'Vdydx ^d^Vdxd/S 
dfidy dx dx dydx dx dx dad^ dx dx 

^dV^ dVd^ dV^ 
da da? dfi dm? dy d^ ' 

d^y ^y 

-j-^ and -T-j being similarly formed, we see that, when the 

three expressions are added together, the terms involving 

dV dV dV 

-J- 9 -ITT, -7- will disappear by the conditions (1), and those 

^•"^^^^^"^ my^x' ^^^ *^^ conditions (2). Hence 
4. Now, let 
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^ r^' d± 

All these expressions satisfy the conditions (1), for a is 
the potential of a homogeneous ellipsoidal shell, of proper 
density, at an external point, and ^ and 7 possess the same 
analytical properties. 

Again, a is independent of v and v\ and is therefore con- 
stant when € is constant. Similarly P is constant when v is 
constant, and 7 is constant when v is constant. Hence a, )8, 
7 satisfy the conditions (2). 

Now 



U) "*■ \dy) ■*■ \dz) 



"(a^ + 6j^6- + €)(C^ + €)iW ■^'Ui^J '^\dz)\' 



2 a 8 

And -2-r' + j3^^-+-2-^- = l• 



with similar expressions for -,- and -p . Hence, squaring 
and adding. 

But from the equations 



a*+e 6* + € c* + € ' a'+v i*-|-v c' + v 

* /3 is a purely imaginary quantity. We may, if we please, write J-lp 
for/3. 
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7 •" It~T~T •" J I ./ "" -■■> 



we deduce 

a? y* ^' _ (® ^ ^) (® — v) (© — V ) 

""a^ + w"" FT© "" c* + ft> "" (o) + a')~(«» + J') (ft) + c') ' 

ft) being any quantity whatever. For this expression is of 
dimensions in m, e, v, v, it vanishes when ft) = e, v, or v\ 
and for those values of ft) only, it becomes infinite when 
ft) = — a*, — J*, or — c*, and for those values of ft) only, and it is 
s= 1 when ft) = X . 

From this, multiplying by a' + a, and then putting 
a = — a', we deduce 

. (e + g") (v + g*) (v + o*) 
(o* - 6'J (a" - c") ' 

a result which will be useful hereafter. 

Again, differentiating with respect to a, and then putting 
e» = e, 

t^ it , "* _ (e — v)(e-v) 

(a" + e/ ■*" (6" + e)* "^ (c' + e)' (e + a") (e + b') (e + C) ' 

. (^\" , (^'\' , fd€\*_. (e + a*){e + V){e + t') 
' ■ \dx) '^\dy) "^ \dz) ~* (e-v) (e-v) ' 



• • W Uy W^/ (e - 1/) (e - v) ' 
(v — u ) (w — e) (e — v) I. aa 






The equation V* 7"= is thus transformed into 
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or {v'-v) [{(e + a')(e + 6') (e + c^}4^JV. 

+ (6 - 1;') [{(v + a*) (v + 5') (v + c*)}* |jj V 

5. A class of integrals of this equation, presenting a close 
analogy to spherical harmonic functions, may be investigated 
in the following manner. Suppose ^ to be a function of e, 
satisfying the equation 

r{(€ + a«)(€ + J')(6 + c")}i^1V=(m€ + r)^, 

m and r being any constants. 

Then, if iff and H^ be the forms which this function 
assumes when v and v are respectively substituted for e, 
the equation V'F= will be satisfied by V^EHH'. 

6. We will first investigate the form of the function 
denoted by E, on the supposition that j& is a rational integral 
function of e of the degree n, represented by 



[ 



We see that 



8' 



= n Kto - 1) (e + a*) (e + 6*) (e + c*) je""* + (n - 2)p,€ 

, (n-2)(n-3) 



+ + P,} 



1.2 



!>,« +-+i', 



11-8 



4] • 



F.H, 8 
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Hence writing 
we see that 



n 



[(» -1)(6» + S/.e* + 3/.e+/J \^+{n- 2) J)., 



.n-8 



{n-2){n-Z) \ 



+ |(e« + 2/;6+/j|e-' + (n-l)f,,e''-^ + ^"-^^^^^-% ,. 



.n-8 



12 Ae^H-- +!>.}• 



Hence, equating coefficients of like poweis of e, we get 

hn - 1) [{n - 2)p, + 3/,} + 1 {(« _ 1) |,^ + 2/j] = n«ip. + r, 
\in - 1) {(^|^)p. + 3 (« - 2)/.p, + 3/.} 

<^(n-l) 



n 



n 



or, as they may be more simply written, 



= rP« 



ELLIPSOIDAL AND SPHEROIDAL HABMONICS. 115 



n 



|(» - 1) f^* -- 2) ^i' + S'lfi j == w^i>i + n 



» {fc^fcl> (n - I) ,. 4. 8 (. - 1)-M + 8 (»-!)/.} 



n (n — 1) 
= 1.2 > P8 + ^^i 



It thus appears that p^ is a rational function of r of the 
first degree, p^ of the second, p^ of the w***, and when the 
letters Pj,p,...p^ have been eliminated, the resulting equa- 
tion for the determination of r will be of the (n + 1)*** degree. 
Each of the letters p^, p^^'Pn ^^^ have one determinate 
value corresponding to each of these values of r; and we 

have seen that wi = n ( n + 5 ) . There will therefore be (n 4- 1) 

values of jE, each of which is a rational integral expression 
of the v^ degree, n being any positive integer. 

7. But there will also be values of E^ of the w* degree, 
of the form 

(e+y)i(e+c')^|e«-'+(n-l)g,e^+ ^"~^^^^;"^^ g,6-'+...+g,.. 

We thus obtain 

{(e + o«)(e+5«)(6 + Oli^ 

= (e + a*)\ (e + V){€ + c') {n-\) le-* + (n - 2) q/^ 

fft-2)(w-3) ,^ , .1 
+ 172 ?«« + - + !?»-.h 

8—2 
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+ (e + o')i (e + b')\ {n - 1) je*^ + (» - 2) y^e""* 

I 

+ (e+a*)* (e+*») (e+c^ (»-l) (n-2) |€--»+(n-3) j^e-* 



(«-2)(w-3) 

j-g q^€ +... + q. 



(n-3)(«-4) 



Hence 



jl (« + *•) (e + c*) + (« + «*) (e + <0 + (« + «*)(« + **)} 

+ i€ + <^ (6+6') (e + (n-l)(»-2)|€«^ + (n-3)jj€-* 

= (me + r) je""* + (n - 1) Jj^ 



(n-lU«-2)_^. ._ 1 



?.e"^ + ...+ J. 



.-. (w-l)(| + n-2) = m, 

(»-l){2o«+|(6« + c')+|(n-2)y.} 

+ (w-l) (»-2){a*+i* + c' + (n-3)jJ = (n-l)my, + r, 
(n - 1) |(^* + oV + a»i«) j.^ + (n - 2) oW qj^ = rj^,. 
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By a similar process to that applied above, we shall find 
that r is determined by an equation of the n*^ degree, and 

that m=(w--l) fn — ^j, and that each of the letters j^, 

9'2* ••?«-! ^s ^ rational function of r. Thus, there will be n 
solutions of the form 

There will also be n solutions of a similar form, in which 

the factors (e 4- c')* (e + a^)*, (e -!• a*)* (e + V)^ are respectively 
involved. Hence, the total number of solutions of the rfi^ 
degree will be 4n + 1. 

8. We may now investigate the number of solutions of 

the degree n + ^ , w being any positive integer. These will 

be of the following forms : three obtained by multipljdng a 
rational integral function of e of the degree n by (e + a*)*, 
(e + i*)*, (e + c*)*, respectively, and one by multiplying a 
rational integi*al function of e of the degree w — 1 by the 
product 

{(€ + a»)(€ + 6*)(€ + c»)}i 

An exactly similar process to that applied above will 
shew us that there will be n-f 1 solutions of each of the 
first three kinds, and n of the fourth. Hence the total number 
of such solutions will be 3 (ti + 1) + n, or 4n + 3, that is 



(« + i) + l. 



To sum up these results, we may say that the total 
number of solutions of the n^ degree is 4n 4- 1, n denoting 
either a positive integer, or a fraction with an odd numerator, 
and denominator 2. 

Similar forms being obtained for H, H\ we may proceed 
to transform the expression EHH' into a function of a?, y, z, 

9. Consider first the case in which 
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Write this under the form 

Then H= (v — ©J (ir — q>j) ... (v— ©J, 

H' = {v^(o^ {v — 6),) ... (i;' — a)„). 

Hence 

EHir=(€-€o^) (v-ft)i) (i;'-a>,)...(6-Q>J(i;-6)„)(v -a)J. 

Now we have shewn (see Art. 4 of the present Chapter) 
that (e — <w J (ir — 6)J (v — a)j) 

= (a>, + a')(«, + 6')(a>,+ c^)f-T^ + T^+:T^-lV 

Each of the factors of EHH' being similarly transformed, 
we see that EHH' is equal to the continued product of all 
expressions of the form 

(«. + a')(» + 5')(a, + c*)(^ + ^+;^^| iV 

> 

the several values of © being the roots of the equation 
ft)'^ + np^ft)^^ + ^^^"^ \ y^ + ... +l)n = 0. 

As this equation has been already shewn to have {n + 1) 
distinct forms, we obtain (ri + 1) distinct solutions of the 
equation V*F=0, each solution being the product of n 
expressions of the form 



a' + ft> 6' + o) c* + (» 

That is, there will be w 4- 1 independent solutions of the 
degree 2w in x, y, z, each involving only even powers of the 
variables. 

10. To complete the investigation of the number of solu- 
tions of the degree 2/1, let us next consider the case in which E 

,=(6+&')i(6+0*{e-'+(n-l)i>,e''^+ ^^\)^^-\ e«-'+...+p^.}. 
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The object here will be to transform the product 

since the other factors will, as already shewn, give rise to the 
product of w — 1 expressions of the form 

^ +-y!_+_g! 1. 

a^ + oD 6* + 0) c^ + a> 

Now, by comparison of the value of a? given in Art. 4, 
we see that 

(€ + 6«)(t; + 5«)(u' + 6»)(6+0(v + c»)(v' + c^ 

= (6«-c«)(i»-a«)(c"-a')(c*-6')yV. 

Hence, we obtain a system of solutions of the form of 
the product of (n — 1) expressions of the form 

a? y" z' ^ 

a+ft> b +(0 c +(0 

multiplied by yz. Of these there will be n, and an equal 
number of solutions in which zx, xy, respectively, take the 
place of yz. 

Thus, there will be 4n + 1 solutions of the degree 2n in 
the variables of which n + 1 are each the product of n 
expressions of the form 

_^_+_y!_+_5!_-i 

a^ + o) 6' + G> c* + ft) ' 

n are each the product of (n— 1) such expressions, multiplied 

byy^, 

W ••• ••• ••• ZXf 

w ••• ••• ••• xy, 

11. We may next proceed to consider the solutions of the 
degree 2n + 1 in the variables x, y, z. 

Consider first the case in which 
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4 

Here the product (e + a')* (u + a*)* (v + a*)* will, as just 
shewn, give rise to a factor x in the product EHH'. 

Hence we obtain a system of solutions each of which is 
the product of n expressions of the form 

a' + CO 6' + © c" + © * 

multiplied by x. Of these there will be 7i4-l, and an equal 
number of solutions in which y, Zj respectively take the 
place of the factor x. 

Lastly, in the case in which 
^ = (6 + a^*(e + J^*(6 + 04J6"-^ + (n-l)p,< 

+ f;2 ^« + — +i^-iN 

we see that in EHH' the product 

(e+a'f (i;+a»)4 (i/'+a^* (6+ J')* (y+b')^ (u' + J')* (e+c')* 

will give rise to a factor an/z. 

Hence we obtain a system of solutions each of which is 
the product of {n — 1) expressions of the form 



.6"-* 



a + CD 6* + © c* + © 

multiplied by a^y^r. Of these there will be n. 

Thus there will be 4n + 3 solutions of the degree 2n + l 
in the variables, of which 

{n + 1) are each the product of n expressions of the form 

9 9 9 

X ti Z 

-J- h Tan ^■ -ft-; 1 multiplied by x, 

(n + 1) are each the product of n such expressions, multiplied 

V** I A^ ••• ••• ••• ••• Zf 



ELLIPSOIDAI^ AND SPHEROIDAL HABMONICS. 121 

n are each the product of {n — 1) such expressions, multi- 
plied by xyz. 

12. Now an expression of the form C . EHH", G being 

any arbitrary constant, is an admissible value of the potential 

a? v' «' 
at any point within the shell -g + '^ 4- -2 = 1. But it is 

not admissible for the space without the shell, since it 
becomes infinite at an infinite distance. The factor which 
becomes infinite is clearly E^ and we have therefore, to 
enquire whether any form, free from this objection, can be 
found for this factor. We shall find that forms exist, bearing 
the same relation to E that zonal harmonics of the second 
kind bear to those of the first. 

Now considering the equation 

which we suppose to be satisfied by putting [7= E, we see 
that, since it is of the second order, it must admit of another 

particular integral. To find this, substitute for U, E jvde, 

we then have 

= [{(€+ a*) (e + 6») (6 + c")}* ^1 E . jvde 

+ (e + a')(e + h')(e + c')^.v 

+ 1 {(« + **) (6 + + (e + c*) (e + a*) + (fi + o») (e + i*)} Ev 
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Now, since by supposition, the equation for the determi- 
nation of Z7 is satisfied by putting U= E, it follows that 

when Elvde is substituted for U, the terms involving \vde 

will cancel each other, and the equation for the determina- 
tion of V will be reduced to 

Tj,dv i^dE 1 f 1 1 1 \ rrl 

•^:j- + i2-t- + s( i + — 15 + 'AE\v=^Oy 

de \ rfe 2VeH-a* € + i* € + c7 J 

Idv 2d^ 1/_1_ _1_ 1 \ 
^^ vde'^ Ede '^2\e + a^^ e + V^ € + (?)' ' 

whence log t; + 2 log ^ + log {(€ + a!) (e + 6*) (e + c*)}* 

= log ^0+ 2 log E^-\- log ahc, 

% and E^ being the values of v and E, corresponding to € = 0. 

Hence t; = v«T^ n; 

••. E Ivde^v^E^abcE \-—- ,, , ^ ,,. , 5771 • 

We may therefore take, as a value of the potential at 
any extemi point, 

F= v,E' abc EHH' f — r . 

For this obviously vanishes when e — oo. It remains so 
to determine v, that this value shall, at the surface of the 
ellipsoid, be equal to the value C, EHB', already assumed 
for an internal point. This gives 

C=v^.E}abc I r. 

Hence, putting v^ . E* . dbc = F,, we see that to the value 
of the potential 

Y.EHH' r — , • 
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for any internal point, corresponds the value 

for any external point. 

13. We proceed to investigate the law of distribution of 
density of attracting matter over the surface of the ellipsoid, 
corresponding to such a distribution of potential. 

Now, generally, if Sn be the thickness of a shell, p its 
volume density, the difference between the normal compo- 
nents of the attraction of the shell on two particles, situated 
close to the shell, on the same normal, one within and the 
other without will be 47r/oS?i. This is the attraction of the 
shell on the outer particle, minus the attraction on the inner 
particle. 

But the normal component of the attraction on the outer 

dV 
particle estimated inwards is —-r- • 

And, if F' denote the potential of the shell on an in- 
ternal particle, the normal component of the attraction on 

dV 

it estimated inwards is ,— . 

an 

Hence ^wpon = -j -j— . 

'^ dn dn 

dV^dVdxdVdydVda 
dn " dx dn dy dn dz dn ' 

dos 
And -J- is the cosine of the inclination of the normal at 
dn 

the point x, y, z to the axis of x, and is therefore generally 
equal to e , , e denoting the perpendicular from the 
centre on the tangent plane to the surface 

^ , y* , ^ _i 
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And we have shewn that 



whence 



or 



2dx _ 1 
a? de " a' + € ' 

^ 2i — • 



da? _ ^ daj 
art ae 



Similarly ^=2e$^. ^=24', 

a/i ae aw ae 

' ' dn \dx de dy de dz de) de 



Similarly -r— = 2e —7- . 

da de 



Now F = K . ^S]ff 



f 

Jo 



(26 



o^»{(6 + a')(e + 6')(6+c')l*' 



de 



And F= V^.EHH' [ — n- ; 

J, IP'{{a* + e){b' + e)i(f+e)}i' 

therefore, generally, 

dV_ Y CTw £^ r de 

de • • deJ,E^ {(a' + e) (b* + e) (c' + e)}* 

But, when the attracted particle is in the immediate 
neighbourhood of the surface, e = 0. Hence, the first line 
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of the value of -— becomes identical with the value of -3—, 

de de 

and we have 

de de ' E, abc' 
Eq denoting the value which E assumes, when € = 0. 

JTTJ' 1 

Hence, ^irphn = 2eV^ -^ — r- . 

But Bn, being the thickness of the shell, is proportional to 
e, and we may therefore write ^ = ^r > ^ being the thick- 
ness of the shell at the extremity of the greatest axis ; 

^ Fq g 1 HE' 
^ 2TrBaabc E^ ' 



• ■ 



and this is proportional to the value of V corresponding to 
any specified value of e, since HH' is the only variable 
factor in either. 

Hence functions of the kind which we are now considering 
possess a property analogous to that of Spherical Harmonics 
quoted at the beginning of this Chapter. On account of 
this property, we propose to call them Ellipsoidal Harmonics, 
and shall distinguish them, when necessary, into surface and 
solid harmonics, in the same manner as spherical harmonics 
are distinguished. They are commonly known as Lamp's 
Functions, having been fully discussed by him in his Legons. 
The equivalent expressions in terms of x, y, z have been con- 
sidered by Green in his Memoir mentioned at the beginning 
of this chapter, and for this reason Professor Cayley in his 
" Memoir on Prepotentials," read before the Royal Society 
on June 10, 1875, calls them " Greenians." 

We may observe that the factor 

1 a 1 



4flr ha abc 
is equal to . , k , axid therefore also to :; kt or -z — r^* 
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Hence, it is equal to 



47r 

-^ {bcBa + caSb + ahSc) 

«> 



or to 



volume of shell ' 
and the potential at any internal point 

/•OO T 

= i volume of shell x EE^ . p I ; r , 

^ ' ^h i?«{(a*+6)(i«H-€)(c'H-€)}* 

and the potential at any external point 

/•OO •• 

= i volume of shell x EE^.p \ r; 

where for p must be substituted its value in terms of i; and v . 

14. We will next prove that if F^, F, be two different 
ellipsoidal harmonics, dS an element of the surface of the 

ellipsoid, UeV^V^dS^O, the integration being extended all 

over the surfaca 

We have generally 

And throughout the space comprised within the limits of 
integration, V» V^ - 0, V» F, = 0. Hence 

//<^.f-''.f)^«-»- 

Now it has been shewn already that Fj, F^ are each of 
the form EHEC^ where jP is a function of e only, H the same 
function of v, R' of u • We may therefore write 

. F.=/.(e)/.M/.(v'). 

and similarly F, =/, (e)/, (w)/, (u). 
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Hence ^i^=^i^,47T. 



...F^-F'^^- 






Now, all over the surface, e = 0. Hence 

/*' fO) f (0) 
Hence, unless -f-j^ — J nJ = 0, which cannot happen 

unless the functions denoted by f^ and /, are identical*, or 
only differ by a numerical factor, we must have 

Now e is proportional to the thickness of the shell at 
any point. Calling this thickness &, we have therefore 



jjSeV,r^d8=0. 



Hence, adding together the results obtained by integrating 
successively over a continuous series of such surfaces, we get 



/// 



V^V^dxdydz=^0; 



F, , Fj now representing solid ellipsoidal harmonics, and the 
integration extending throughout the whole space comprised 
within the ellipsoid. 

* This may be shewn more rigorously by integrating through the 
space bounded by two confocal ellipsoids, defined by the values X and /a of c. 
We then get, as in the text, 

Now the factor within { } eannot yanish for aU values of X and fi, unless the 
functions devoted by fi and /, be identical, or only differ by a numerical 
factor. 
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15. It will be well to transform the expression 



II' 



to its equivalent, in terms of v, v\ 

For this purpose we observe that if ds, ds' be elements of 
the two lines of curvature through any point of the ellipsoid, 
dS= dads'. 

Now, 
da^ is the value of ds^ + dy^ + dz^ when e and v are constant, 
ds^ ... • • • ... 6 and v 

and <r»- (^ + «')(v+«')('^' + «') . 

*°^ '^~ {b* - a^ {fi* - a*) 

therefore if e and v do not vary, 

2dx dv 

••. dw = Ti-^ du. 

2a +v 

Similarly c?y = |^^rfv, d^ = -^dv; 

1 f >^ V* «* 1 



Again, differentiating with respect to © the expression 

stined for -n \- tt- — + -^ 1> wo get 

a^'H-o) 6'+o) c*H-a) ^ 



aj" , ^ , 2* . (v — ©) (u — ©) 



{y — co) (e — o) (c — ©) (v — ©) (v' — ©) 
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therefore, putting a — v, 



X* 



f . _l {V'-V ){€-V) . 



(o* + v)» "'" (6' '+ u)' "'" (c» + 1/)* " (<i» + u) (J* + y) (cP + v) ' 

,,_1 {v —v)(e — v) J , 
•*• *" ~ 4 (a» + v) (J*+ v) (c' + v) '^''^' 

A similar expression holding for ds'* we get 

jo.__JL (v-v)*(e-v)(6-v) , ,, „ 

" 16 (o'+u) (y +w) (c'+u) (a»+ v') (6'+v') (c'+v') '*" ' 

A • 1 _ a^ . y* . ^^ 

Again, g.-(ji.^.^). + (j. + e). + (c« + e)» 

_ (e-v')(e-„) 
~(a'+e)(J' + 6)(c' + €)' 

writing e for q> in the expression above ; 

It has been shewn that, integrating all over the surface, 
the limits of v are - c' and —J', those of v , — i* and — a*. 

Hence, F^, F,, denoting two different ellipsoidal har- 
monics 

]-i^]-a^ {(a»+i;) (6«+v) (c'+i;) (a»+ v ) (6*+u ) (c'+i;')}* 

The value of the expression \\\ V^dxdydz, or its equiva- 
lent 

, f-^ /•-»• V^{v'^v)dvdi/ 

^ j-6*j-a«{(a»+i;)(i»+i;)(c»+v)(a*+i/')(t'+i'')(c*+i;'j}*' 

in any particular case, is most conveniently obtained by 
expressing F as a function of x, y, z, 

F.H. 9 
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16. Before proceeding further with the discussion of ellip- 
soidal hannonics in general, we will consider the special case 
in which the ellipsoid is one of revolution. We must enquire 
what modification this will introduce in the quantities which 
we have denoted by a, A 7> viz. 



« = f- 

^c la 



df 



B=r ^ 

and in the differential equation 



We will first suppose the axis of revolution to be the 




J, e'-a*+b* (a*-b*)i ^ ,-(o«-6«)4' 

~ = 2f" «^^ _ 1 , (a'-h')i-<o 
^ Jo 0'-a*+ V - ia^^h^h '*'S (a. _ j.)* + „ ' 

To transform ^, we mast proceed as follows. 

Put -i/rsr-c'cosV-J'sinV, i; = -c*cos*^-6»sin'<^, 
we then get generally 

6« + ^=(6«_c')cosV, c' + ^=(c'-6»)sinV; 

dy^ = 2 (c'* — 6*) cos «r sin «r d!«r ; 

• B- ^ r ^"^ ^ 2V^ 
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Hence, 



i..-i(,._.vii|, 



Also, e = t}' — a', v' = (o*-a', v = — V, and our differential 
equation becomes 

+ (i,' -a^ + V) ko* -a*+ 5') £r F 
_(<„«_,»)(a«-J«)^=0, 

-(a'-6')(,«-a,')^=0. 

This equation may be satisi&ed in the following ways. 

First, in a manner altogether independent of <^, by sup- 
posing V to be the product of a function of 97 and the same 
function of o), this function, which we will for the present 
denote by/ (97) or/(cD), being determined by the equation 



or 



or 






9-2 
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Secondly, by supposing ^ -^ a constant multiple of F, 
= — cr*F, suppose. 

Our equation may then be written 

(o,»-a« + 6')|(,'-a' + 6')|}V 

_(,'_a' + J«)|(o,»_«'+6«)^JV 

- ff* (a' - V) {(«■ - a» + 6*) - {if - o» + i*)} F= 0, 

which may be satisfied by supposing the factor of V inde- 
pendent of ^ to be of the form F{^ ^(«)» where 

|(,«_ o» + V) ^\'f(v) - o* (a* - 6') F{7,) = m iv'-a*+b') F{v), 

y-o' + ¥) ^}V(«) -«7» {a*-b*) F{o>)=m{a*-a*+h') F{ti>). 

The factor involving ^ wiU be of the form 

A cos a<j)+JB sin o-^. 

Now, returning to the equation 
d 



{ 



drj 



■(,•- a' + 5») II /(,) = «./(,), 



we see that, supposing the index of the highest power of rj 
involved iDifiv) to be i, we must have m = i(i+ 1). 

Now, it will be observed that 97 may have any value 
however great, but that ©*, which is equal to a* + v, must 
lie between a* — 6* and 0. Hence, putting ©' = (a* — i') /a', 
where fi* must lie between and 1, we get 
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Hence this equation is satisfied by /{(a*— b*)^ fi} = CP„ 
G being a constant ; and supposing u = 1 we obtain the 
following series of values for/ (<»), 

t = o, /w = i, 

t = l, /(a>)= " 



i=2, /(«) 



2 (a" - b") 



^\ ' 



Exactly similar expressions may be obtained for /(rj), and 
these, when the attraction of ellipsoids is considered, will 
apply to all points within the ellipsoid. But they will be 
inadmissible for external points, since 17 is susceptible of in- 
definite increase. 

The form of integral to be adopted in this case will be 
obtained by taking the other solution of the diflferential 
equation for the determination of /{v), i»e. the zonal har- 
monic of the second kind, which is of the form Q. \ rl- , 

where 

d9 






.(a — V)' 

Or, putting 17' = (a* — J*) i/', ^ = (a* — i") V, we may write 

d\ 



Q,{v)^PAv)[ 



^*Wi(^-i) 



17. We may now consider what is the raeaninof of the 
quantities denoted by rj and c». They are the values of 9r 
which satisfy the equation 
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and are therefore the semi-axes of revolution of the surfaces 
confocal with the given ellipsoid, whicL pass through the 
point X, y, z. One of these surfaces is an ellipsoid, and 
its semi-axis is 17. The other is an hyperboloid of two sheets 
whose semi-axis is co. 

Now, if be the eccentric angle of the point x, y, z, 
measured from the axis of revolution, we shall have 

a? = 17' cos* 0. 

But also, since 17*, a>\ are the two values of y which 
satisfy the equation of the surface, 

Hence ©* = (a* — 6") cos' 0, 

and we have already put 

whence the quantity which we have already denoted by fjk 
is found to be the cosine of the eccentric angle of the point 
X, y, z considered with reference to the ellipsoid confocal 
with the given one, passing through the point a?, y, z. We 
have thus a method of completely representing the potential 
of an ellipsoid of revolution for any distribution of density 
symmetrical about its axis, by means of the axis of revo- 
lution of the confocal ellipsoid passing through the point 
at which the potential is required, and the eccentric angle 
of the point with reference to the confocal ellipsoid. For 
any such distribution can be expressed, precisely as in the 
case of a sphere, by a series of zonal harmonic functions of 
the eccentric angle. 

18. When the distribution is not symmetrical, we must 
have recourse to the fonn of solution which involves the feictor 
A cos cr<^ + J5 sin cr^. It will be seen that, supposing F to 
represent a function of the degree i, and putting m = i (4+I), 
the equation which determines F(oi) is of exactly the same 
form as that for a tesseral spherical harmonia For -^^(i;), if 
the point be within the ellipsoid, we adopt the same form, 
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if without, representing the tesseral spherical harmonic by 
r/') I — 5— jl , or yw (v), we adopt the form 



2'/"('')£ 



d\ 



5FWl'(V-i)* 



19. It maybe interesting to trace the connexion of sphe- 
rical harmonics with the functions just considered. This may 
be effected by putting 6' = a'. We see then that rj will become 
equal to the radius of the concentric sphere passing through 
the point, and rf — ci?-\-V will become equal to rj^. Hence 
the equation for the determination o{f{rj) will become 

d / ^ d 



|(.-|)/(.) = «(»4- !)/(.). 



which is satisfied by putting f{rj) =rf, or rf^**^K The former 
solution is adapted to the case of an internal, the latter to 
that of an external point. 

With regard to f(co), it will be seen that the confocal 
hyperboloid becomes a cone, and therefore <a becomes inde- 
finitely small. But fi, which is equal to 1, remains 

(a* — 6')* 

OS 

finite, being in fact equal to - or cos 0. Hence /(/i) becomes 
the zonal spherical harmonic. 

Again, the tesseral equations, for the determination of 
F{r}),F{(o), become 

which are satisfied by Ffjj) =rf or 97"^'*^*. 

And, writing for <o\ (a*- J') fi\ we have, putting F((o) =x(/^)* 

(m'-i)^}'xW+<^x(a*)=«(»'+i)G^*-i)xW. 

which gives % (/i) = 3'/*'^ (/a). 
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. 20. We will next consider the case in which the axis of 
revolution is the least axis of the ellipsoid, which is equi- 
valent to supposing a*=b\ To transform a and )8, put 
c* + -^ = ^*, c^ + € = r}\ c' + 1; = o)', we thus obtain 

„ _ 2 r _^?_ -_ A_ ten- («*-^ 

To transform 7, we must proceed as follows : 

Put -^ss — a'sin'«r — y cos'o-, v' = — a' sin* ^ — 5' cos' ^, 
we then get, generally, 

a' + i|r=(aV-i*)cos»i3-, 6» + '^fr= -(a"-6») sin'isr, 

c'+'^=c'— a'sin'^— J*cos*^, e?^=— 2 (a*— 6") sinw costar rfw. 

Hence 

'^ i« (a* sinV + 6" cosV - O* (a" - c*)* 
Hence, |- = -|(a«-c' + ,«) ^. 

_ = __(« _c +««)-. 
also, 6 = 97* — c*, 



3 s 
V=0) — c, 

' a 

V = — a , 



and our differential equation becomes 

dV 



-(a»-c' + i7')|(a'-c' + a>«)4}V 
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We will first consider how this equation may be satisfied 
by values of F independent of (j). 

We may then suppose Fto be the product of a function 
of V, and the same function of o), this function, which we will 
suppose to be of the degree i, being determined by the 
equation 






|(a'-c'+a,»)|^|/(«,)=i(« + l)/(«). 



On comparing this with the ordinary differential equa- 
tion for a zonal harmonic, it will be seen that, on account 
of a' being greater than c", the signs of the several terms in 
the series for f(rj) will be all the same, instead of being 
alternately positive and negative. We shall thus have 

* = 1. fiv) = " 



»=2, /(,) = 



(a'-c')i' 
STi' + a^-c? 



2 (ft» - c') 



,«\ » 



.-3. /(,).»?^±i(^, 

2 (a* - c')' 

. _ B5v* + 30 (a* - O 1?' + 3 {a* - <?)' . 
»-*, /W 8 (a' - O' ' 



and generally 



•^("^ = (^4 274:^ I' ("'+"*- '''^' 

f V 
We will denote the general value of /(ly) by jp, \ — ^ — — rj- , * 

or, writing 17 = (a' — <f)^v, by Pi {p). 
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For external points, we must adopt for / (17) a function 
which we will represent by qA , , x r > or 2,(1/), which will 
be equal to 



* l(a»-c*)ij 



(^+a»-c^ 



It is clear that /(w) may be expressed in exactly the 
same way. But it will be remembered that 47* and to* are 
the two values of y which satisfy the equation 

^ + 2^+^ = 1. 



Hence 17, as before, is the semi-axis of revolution of the 
confocal ellipsoid passing through the point (a?, y, z). But 
97*0)* = — (a* — c') «*, an essentially negative quantity, since 
a* is greater than c". Hence o)* is essentially negative. Now, 
if ^ be the eccentric angle of the point {x, y, z) measured 
from the axis of revolution, we have «' = 17' cos'^. Hence 

^V = -(a»-c»)i7'cos"^, 

and therefore ©' = — (a' — c') cos' 

= — {a^ — c*) fi^, suppose. 

Hence the equation for the determination of /(«) assumes 
the form 

the ordinary equation for a zonal spherical harmonic. Hence 
we may write 

/(«)=p.H 

fi being the cosine of the eccentric angle of the point x, y, z, 
considered with reference to the confocal eUipsoid passing 
through it. 
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21. We have thus discussed the form of the potential, 
corresponding to a distribution of attracting matter, sym- 
metrical about the axis. When the distribution is not 
symmetrical, but involves <^ in the form A cos a'<f> + B sin <r^, 
we replace, as before, P, (/i) by T/''^ (/i), and Pi(fi) by a 
function t^^'^^(y) determined by the equation, 

and a, (v) by ^/'> (v) I , , ■„ — r . 

22. As an application of these formulae, consider the fol- 
lowing question. 

Attracting matter is distributed over the shell whose 

2 7/* _L «* 

surface is represented by the equation -7+—r2~ — ^f so 

that its volume density at any point is P< (m), /a being the 
cosine of the eccentric angle, measured from the axis of 
revolution ; required to determine the potential at any 
point, external or internal. 

The potential at any internal point will be of the form 

CPMPiiy) (1), 

and at an external point, of the form 

C'PMQM (2), 

where (a* — V)^ v = the semi-axis of the figure of the con- 
focal ellipsoid of revolution passing through the point (/i, v). 

Now the expressions (1) and (2) must be equal at the 

surface of the ellipsoid, where v = --- — ---. . 

(a* - i')* 

Hence 
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But generally 

Hence 

dK. 



'^'V-i')*}" '^V-6'j4jL 



»_i>.(\)l'(X»-l)' 



d\ 



We may therefore, putting C = AP^ \—^ ^J- , write 



and we thus express the potentials as follows : 

^-P« (m) Pi (y) Qt I ,\a^ r at an internal point, 

AP^ (ji) Q^ (p) P^ \ — — ^V at an external point. 

Or, substituting for Q^ its value in terms of P,, 



^^ ^' ^"^ ^' {(7r^,i} /Li-,^^' (V - 1) 

at an internal point, 

at an external point. 

Now, to determine -4, we have, Sa being the thickness 
of the shell at the extremity of the axis of revolution, 
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^ 4!tr Ba , rf \ drj drj Jri=a 

1 



1 a A p, V 1 
Si- Sa a* - 6'- ' ^^ "S* 7 

Ts ~ 1 



^'V-w W"%"0 



a'^b' 



■"Z1T2£:-4P, (/i). 



47r 6*&t 
Hence, if p = P, (/a), we obtain 



^ = 47r— =47rJ8J. 
a 



And we thus obtain 



dX 



(a«-6a)4 



i>.wr(x.»-i) 



K = 4Tr6S6 P, (^) P, {,-r^— i) P. M r ^==^ 

= 4^6SJP.(M)Q,(.)P,{^^f^}. 

If the shell be represented by the equation 

it may be shewn in a similar manner that we shall have 
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Fj = 4firada P, (/*) p, {v) q, j _^^^ | , 

r, = 4i-irada P^ (/i) q,{v) p^ | ^ X . 

23. We may apply this result to the discussion of the 
following problem. 

If the potential of a shell in the form of an ellipsoid of 
revolution about the greatest point he inversely proportional 
to the distance from one focm, find the potential at any 
internal point, and the density. 

If the potential at P be inversely proportional to the 
distance from one focus S, and H be the other focus, we have, 

HF+SP = 2rj, HP-SP=^2a>, 

.'. SP = 'n'-<o. 

Hence if M be the mass of the shell, F, the potential at 
any external point, ' 



M 



M 



(a*-6»)i 



S (2i+l) P. 0*) <?.(,;). 



Now, by what has just been seen, the internal potential, 
corresponding to P, (ji) Qt (y), is 

Pi W Pi (") — j — - — 1^ 

Hence, if F, be the potential at any internal point. 
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And the volume density corresponding to P, (/it) ^^ (y) is 

4^686 pi — £_].* 

Hence the density corresponding to the present distri- 
bution is 



P = ^^^-^ — S (2t + 1) 



*W-6»)4J 



.(a*-6')*J 
If F, had varied inversely as HP, we should have had 

and our results would have been obtained from the foregoing 
by changing the sign of ©, and therefore of /a. 

24. Now, by adding these results together, we obtain 
the distributions of density, and internal potential, corre- 
sponding to 

* 17 — ft) 17 + a) 1^—0)* 

or, in geometrical language, 

^^" SF'^HP 8F.BF ' 

= M multipHe<l by the axis of revolution of the confocal 
ellipsoid, and divided by the square on the conjugate semi- 
diameter. We may express this by saying that the potential 
at any point on the ellipsoid is inversely proportional to the 
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square on the conjugate semi-diameter, or directly as the 
square on the perpendicular on the tangent plane. 

Corresponding to this, we shall have, writing 2k for *, 
since only even values of i will be retained, 



f.' 



2M 



Q 



ik 



V - *')') 



'^■-r^3yit^(*^-^') „ ,-r /^j f.Wi-.W. 



|{a' - m 



h being 0, or any positive integer. 

Again, subtracting these results we get 



* 1^ — 0) 17 + 0) i;*— ft)' * 



= ilf multiplied by the distance from the equatoreal plane, 
and divided by the square on the conjugate semi- diameter. 

This gives, writing 2A; + 1 for i, 



<?. 



w - m 



25. In attempting to discuss the problem analogous to 
this for an ellipsoid of revolution about its least axis, we see 
that since its foci are imaginary, the first problem would re- 
present no real distribution. But if we suppose the external 
potential to be the sum or difference of two expressions, each 
inversely proportional to the distance from one focus, we 



ELLIPSOIDAL AND SPHEBOIDAL HABMONICS. 145 

» 

obtain a real distribution of potential — in the first case 
inversely proportional, to the square on the conjugate 
semi-diameter, in the latter varying as the quotient of the 
distance from the equatoreal plane by the square on the 
conjugate semi-diameter. 

It will be found, by a process exactly similar to that just 
adopted, that the distributions of internal potential, ia»nd 
density, respectively corresponding to these will be : 



In the first case 



2M 



'Ai^^m 



Jc being 0, or any positive integer. 
In the second case 



P= . ,r\i . S(4A; + 3) 



Z: being 0, or any positive integer. ^ 

26. We may now resume the considei:ation of the ellip- 
soid with three unequal axes, and may shew how, when the 
potential at every point of the surface of an ellipsoidal shell 
is known, the functions which we are considering may be 
employed to determine its value at any internal or external 
point. We will begin by considering some special cases, 

F. H. IQ 
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by which the general principles of the method may be made 
more intelligible, 

27, First, suppose that the potential at every point of 
the surface of the ellipsoid is proportional to a? = -^— suppose. 

In this case, since x when substituted for V, satisfies the 

equation V'^= 0, we see that V^- will also be the potential 

at any internal poii^t. But this value will not be admissible 
at external points, since a becomes infinite at an infinite 
distance. 

Now, transforming to elliptic co-ordinates 



(C 



'}' 



And the expression 
F, f (e+g') (v+o') (v'+o') | I r d-^ 

. r . d'^ 

■ io (f + a") {(t + a") (t + V) (f + c')]* 

satisfies, as has already been seen, the equation v'F"=0, is 

equal to Fj,— at the surface of the ellipsoid, and vanishes 

at an infinite distance. This is therefore the value of the 
potential at any external point. It may of course be written 

v^r d^ 



a J, 



dy^ 



4 



28. Next, suppose that the potential at every point of 
the surface is proportional to yz^VA-^ suppose. In this 
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case, as in the last, we see that, since yz when substituted 
for V, satisfies the equation ^^V=^0, the potential at any 

internal poitit will be Vojrl while, substituting for y, z their 

values in terms of elliptic co-ordinates we obtain for the 
potential at any external point 

y,y^ r df 

dyfr 



^. 



29. We will next consider the case in which the po- 

tential, at every point of the surface, varies as a?' = Fo -5 

suppose. This case materially differs from the two just con- 
sidered, for since a? does not, when substituted for V, satisfy 
the equation v*F= 0, the potential at internal points cannot 
in general be proportional to x\ We have therefore first to 
investigate a function of x, y, z, or of 6, v, v which shall 
satisfy the equation y^F==0, shall not become infinite within 
the surface of the ellipsoid, and shall be equal to oo^ on its 
surface. 

Now we know that, generally 

- (a' + o)) (6* + ©) (c» + ft)) = (€ - 0)) (1/ -0)) {v - o)). 

And, if ^^, ^jj be the two values of <o which satisfy the 
equation 

(6«+a)) (c'+o)) + (c'+o)) (a'+o)) + (a'+(») (6'+ft))=0...(l), 

we see that 

V»(e-^,)(v-^0(u'-^.)=O. 

and ^«(e_^^(u_^J(u'-^J = 0. 

And, by properly determining the coefficients A^, A^, A^, 
it is possible to make 

A,+A,(e-e;)(v-0;)(y'-0;)+A,{e-0;){v-e,){v'-0;)...{2) 

= ^ when fcV + cVw* + a'6V - o*6V = 0. 

IQ— ^ 
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Hence, the expression (2) when A,, A^, A^ are properly 
determined will satisfy all the necessary conditions for an 
internal potential, and will therefore be the potential for 
every internal point. 

Now, we have in general 

and, over the surface 

Pc V + c*ay + a»6V - a'6V = 0. 

Hence, S^ being any quantity whatever, we have, all over 
the surface, 

and therefore, putting ^ = — a', 

(a«-6') (a»-c»)x'=^^±|j (e-^O (v-f,) («'-^J 

Hence, the right-hand member of this equation possesses 
all the necessary properties of an internal potential. It 
satisfies the general differential equation of the second order, 
does not become infinite within the shell, and is proportional 
to ic* all over the surface. 

We observe, by equation (1), that 
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identically, and therefore, writing — a* for «t, 

(a* - h*) (a' - c*) = 3 (a* + 6^) (a* + ^J. 

Hence, over the surface of the shell. 



scr = 



a* 



and we therefore have, for the internal potential. 



This is not admissible for external points, as it becomes 
infinite at an infinite distance. We must therefore substi- 
tute for the factor e — ^^ 

Ce-e) r ■ ^ ■ I 

■ J. (^ - dj" {(t + a') {-^ + V) (ir + c")]* 

with a siokilar substitution for e — dg, thus giving, for the 
external potential, 

Tr-n p-g.)(»-g.)(v-g.) r d± 

; ^f i±_ I 

■ Jo (t - ex {(^ + a') (f + 6') (f + c') jl 
~eJ^F^J{^+eS J. (^-^/ {(t+a') (^+6^) (f +c*)l* 
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The distribution of density over the surface, corresfpopd- 
lug to this distribution of potential, may be investigated by 
means of the formula ' 

^ 27rda\ de de /,so* 
or its equivalent in Art. 13 of this Chapter. We thus find that 
1 g F, r (v-g,)(v--g.) 






d^jr 



+i^r *ti^ 1 

30. The investigation just given, of the potential at an 
external point of a distribution of matter giving rise to a 
potential proportional to a? all over the surface, has an in- 
teresting practical application. For the Earth may be re- 
garded as an ellipsoid of equilibrium (not necessarily with 
two of its axes equal) under the action of the mutual gravi- 
tation of its parts and of the centrifugal force. If, then, 
y denote the potential of the Earth at any point on or with- 
out its surface, and il the angular velocity of the Earth's 
rotation, we have, as the equation of its surface, regarded as 
a surface of equal pressure, 

/. F 4- ^ fl* (a!^ + ^) = a constant, 11 suppose. 

Hence, if a, h, c denote the semi-axes of the Earth, we 
have, for the deterinination of F, the following conditions ; ' 
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S-f-^"-»-" •«. 

y=0 at an infinite distance (2)| 

F=n-|il»(a^4-y') when • 

-ii + ^ + -i = l....^ (3). 

or b c 

The term 11 will, as we know, give rise to an external 
potential represented by 

nf / ^^ ■ r ^^ 

1 1 . . . ' 

The two terms —^ ilV, — ^iiy, will give rise to terms 

which may be deduced from the value of F^ just given by 

1 1 . 

successively writing for Y^ — ^ilV, and —^11*6', and (in 

the latter case) putting 6' for a' throughout, . We thu« get 
(g-g.)(v-g.)(v-g.) f d^ ■ 

^f ^^L 
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31. Any rational integral fonction V of x, y, z^ -which 

satisfies the equation jfV^ 0, can be expressed in a series 

of Ellipsoidal Harmonics of the degrees 0, 1, 2...» in x^ y, z. 

For if F be <^ the degree i, the number of terms in V will 

, (♦4-l)(*+2)(t + 3) ^ ^, ... .^ ^ . 

be ^ ^ • Now the condition y*F=0 is 

equivalent to the condition that a certain function of x, y, z 
of the degree i— 2, vanishes identically^ and this gives rise 

to ^^ ^-^ conditions. Hence the number of inde- 

D 

pendent constants in F is 

(t + 1) (t-f 2) (t> 3) (t - 1) t (♦ -f 1) 
6 6 ' 

or {i + 1)\ And the number of ellipsoidal harmonics of the 

1 3 % 

degrees 0, 1, 2...» in w, y, z or of the degrees 0, ^, 1, ^-..^ 

in 6, V, V y is, as shewn in Arts. 6 to 10 of this Chapter, 

1 + 3+5+.,. + 2i + l, 

or (i + 1)*. Hence all the necessary conditions can be satis- 
fied. 

32. Again, suppose that attracting matter is distributed 
over the surface of an ellipsoidal shell according to a law of 
density expressed by any rational integral function of the 
co-ordinates. Let the dimensions of the highest term in this 
)Bxpression be f, then by multiplying every term, except those 
of the dimensions i and t — 1 by a suitable power of 

we shall express the density by the sum of two rational inte- 
gral functions of a;, y, z of the degrees i^ i — 1, respectively. 
The number of terms in these will be 

(i±lHi±2)^.l(i^±)^(,+l). 
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And any ellipsoidal surface harmonic of the degree i, i *- 2... 
in a;, t/, z, may, by suitably introducing the factor 

be expressed as a homogeneous function of a?, y, z of the 
degree i ; also any such harmonics of the degree t— 1, t— 3... 
in a?, y, z may be similarly expressed as a homogeneous 
function of aj, y, z of the degree i — 1, And the total number 
of these expressions will, as just shewn, be (t4- 1)*, hence by 
assigning to them suitable coefficients, any distribution of 
density according to a rational integral function of x, y, z 
Tcs2,y be expressed by a series of surface ellipsoidal harmonics, 
and .the potential at any internal or external point by the 
corresponding series of solid ellipsoidal harmonics. 

33. Since any function of the co-ordinates of a point on 
the surface of a sphere may be expressed by means of a series 
of surface spherical harmonics, we may anticipate that any 
function of the elliptic co-ordinates v, u' may be expressed by 
a series of surface ellipsoidal harmonics. No general proof, 
however, appears yet to have been given of this proposition. 
But, assuming such a development to be possible at all, it 
may be shewn, by the aid of the proposition proved in 
Art. 15 of this Chapter, that it is possible in only one way, 
in exactly the same way as the corresponding proposition 
for a spherical surface is proved in Chap. IV. Art. 11. 

The development may then be effected as follows. De- 
noting the several surface harmonics of the degree i in a?, y, «, 

or \ in V, vy by the symbols V,^\ F/^), ...F/«+^), and by 

F{v, v) the expression to be developed, assume 

j'(v,vO= c;,F,+ a/i>F/i)+ a,<2) f,(2)+ (7wf/8)+... 

Then multiplying by cF/''^ and integrating all over the 
surface, we have 

jeF {v, v') F/-> dS = Cl^^je (F/'>)« dS. 
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The values 



o{jeF{v, i/^F; 



<') dS, and 



oijeiVy^f 



(ZiS^must 



be ascertained by introducing tbe rectangular co-ordinates 
X, y, z, or in any other w^y which may be suitable for the 
particular case. The coefficients denoted by C are thus 
determined^ and the development e£fected. 



EXAMPLES. 



1 • ■ * * 1 # 



1. Prove that (sin <?)« = A p^_1|p, + ^P^. . 

Why cannot (sin $)* be expanded in a finite series of spherical 
harmonics? 

1 1 1 1 + sin ^ 

2. Prove that 1 + o ^i+ 3 ^. + 4 A+ - =log — T 



8in ^ 



3. Establish the equations 



4. If ft = cos tf, prove that 



P.(/^)=l-»(*H-l)8in«? + ... + (-ir(|i^^(sin«^) 



+ ... 

:5 " ' {\m)- \i-m \ z/ 

and also that 

i'.Ot)=(-iy+(-i)'+>»(i+i)coB'|+ 



• •• 



V ' ( w)' tr-m\ 2/ 



(|m)'[*- 

5, Prove that, if a be greater than c, and i any odd 
integer greater than m, 

f * / . « • — -* 1 2 ft + m c« 



6. Prove that f rp\ dfi = f (i + 1), 
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.^p It + m (a1)"+* 
7. Prove that, when tt = * 1, -7-5 = . \^/ . 



ci/A" 



8. Prove that 



IP P 



P P P 



«•! 



S 



••• ••• ••• ••• 

P P P 



is a numerical multiple of 



9. Prove the following equation, giving any Laplace's co- 
efficient in terms of the preceding one : 



P^. = pP. + ^f,P.^P'^0, 



where Cp= jiix + JT^^^Jl -fi" cos (co - w') and C is zero if n be 
even, and 



(-i)« 



2'*' mn + l)\ 



-, , if 71 be odd. 



10. If i, j, k be three positive integers whose sum is even, 
prove .that 



/ P^PA^i^ 



2.4...({+y + ;fc) 1 



1.3 ... (i + j* + J5;-l) t+J+A;+l ' 

Hence deduce the expansion of P^Pj in a series of zonal 
harmonics. 

11. Express afy + f^ + yz+y + z as a sum of spherical 
harmonics. 

12. Find all the independent symmetrical complete harmonica 
of the third degree and of the fifth negative degree. 

13.' .Matter is distributed in an indefinitely thin stratum over 
the surface of a sphere whose radius is unity, in such a manner 
that the quantity of matter laid on an element .{SS) of the surface 

is 8/S'(l+aa)+6y + c«+/a;' + ^y*-tfo'), 



. y 
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where x, y, z are rectangular co-ordinates of the element SS re- 
fen*ed to the centre as origin, and a, b, c, /, g, h are constants, 
j^ind the value of the potential at any point, whether internal 
or external. 

14. If the radius of a sphere be r, and its law of denaitj be 

p = ax+b7/ + cZy where the origin is at the centre, prove that its 

47rr* 
potential at an external point (^, rj, C) is ,- „3 (a^+ 617 + c{) where 

It is the distance of (f, rj, ^) from the origin. 

15. Let a spherical portion of an infinite quiescent liquid be 
separated from the liquid round it by an infinitely thin flexible 
membrane, and let this membrane be suddenly set in motion, 
every part of it in the direction of the radius and with velocity 
equal to S^, a harmonic function of position on the surface. Find 
the velocity produced at any external or internal point of the 
liquid. State the corresponding proposition in the theory of 
Attraction. 

16. Two circular rings of fine wire, whose masses are Jf and 
Jtr, and radii a and a', are placed with their centres at distances 
6, 6', from the origin. The lines joining the origin with the 
centi'es are perpendicular to the planes of the rings, and are in- 
clined to one another at an angle 0. Shew that ^e potential of 
the one ring on the other is 



j/jr2:::r(^^.^.«.). 



where B -h'- "(*'~^> 6*-«« + n(n-l)(n-2)(n-5) ^ . 
Where £,-b j-^b a+ 2.2.4.4 * " " - 

and B^^ and Q^ are the same functions of V and a' and of cos 
an d sin respectively, and c is the greater of the two quantities 

Va* + 6» and J^^TV\ 

17. A uniform circular wire, of radius a, charged with 
electricity of line-density e, surrounds an uninsulated concentric 
spherical conductor of radius c; prove that the electrical density 
at any point of the surface of the conductor is 

the pole of the plane of the wire being the pole of the harmonics. 
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18. Of two Bpterical conductors, one entirely surrounds tliei 
other. The inner has a given potential, the outer is at the 
potential zero. The distance between their centres being so 
small that its square may be neglected, shew how to find the 
potential at any point between the spheres, 

19. If the equation of the bounding surface of a homo- 
geneous spheroid of ellipticity c be of the form 



r=a(i-l^,), 



prove that the potential at any external point will be 

M C-A 



r r* 



p,» 



where C and A are the equatoreal and polar moments of inertia 
Of the body. 

Hence prove that V will have the same value if the spheroid 
be heterogeneous, the surfaces of equal density differing from 
spheres by a harmonic of the second order. 

20. The equation E = a(l+ay) is that of the bounding 
surface of a homogeneous body, density unity, differing slightly 
in form and magnitude from a sphere of radius a; a is a 
small quantity the powers of which above the second, may be 
tieglected; and y ia a function of two co-ordinate angles, such 
that 

where T^y T ,,. Z^^ Z^,,. are Laplace's 'functions. Prove that 
the potential of the body's attraction on an external particle, 
the distance of which from the origin of co-ordinates is r, is 
given by the equation 



r 
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21. If Jf be the mass of a uniform hemispherical shell of 
radius c, prove that its potential^ at any point distant r from 
the centre, will be 

M \M(\ -, 1 „»^ 
2c 2 c»V2' 'i^A^c^ 



2.4.6 



^*? 2,4.6. 8^^c«^-7' 



or 



2r 2 \2 V 2.4 V 



, 3 c' 3.5 c' \ 

+ 27476 '^V 2. 4.6.8 '^V'^ "V' 

according as r is less or greater than c; the vertex of the hemi- 
sphere being at the point at which fi = 1« 

I 

22. A solid is bounded by the plane of ajy, and extends to 
infinity in all directions on the positive side of that plane. 
Every point within the circle a* + y' = a*, « = is maintained at 
the uniform temperature unity, and every point of the plane xy 
without this circle at the uniform temperature 0. Prove that, 
when the temperature of the solid has become permanent, its 
value at a point distant t from the origin, and the line joining 
which to the origin is inclined at an angle -Q to the axis of z will 
be 

^»"^»a**"2^*a»""2~4^'a*"*'- 

1.3...(2f-l) 1^ 
^^"^^ 2.4...2i -^«^i a^» •*" ••' 

if r < a, and 



Ip ^ 1.3 



0^ 1.3...(2t-l) a« 



2 ^r» 2.4 »r* ^ ' 2,4...2i 

if r > a. 

23. Prove that the potential of a circular ring of radius c, 
whose density at any point is cos mi^, c^ being the distance of the 
point measured along the ring from some fixed point, ia 
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2a-co8m^(8intf)*ir— ^— ^ — "• 



2m dftr f^*' 
2.4.6...(2m + 2) .c^/x" r' 



+ 



'*"2.4.6...2(w + ;fc diiT r^+»'*"-j' 

wliere r is greater thcui c. If r be less than c, r and c must be 
interchanged. 

24. A solid is bounded by two confocal ellipsoidal surfaces, and 
its density at any point P varies as the square on the perpendicular 
from the centre on the tangent plane to the confocal ellipsoid 
passing through P. Prove that the resultant iittraction of such 
a solid on any point external to it or forming a part of its mass 
is in the direction of the normal to the confocal ellipsoid passing 
through that point, and that the solid exercises no attraction on a 
point within its inner surface. 
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